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TECHNICAL NOTE NO. 1832 

SMALL BENDING AND STRETCHING OF SANDWICH-TYPE SHELLS 
Ity" Eric Reissner 

SUMMARY 


A theory has been developed for small bending and, stretching of 
sandwich-type shells. This theory is an extension of the known theory 
of homogeneous thin elastic shells. It was found that two effects are 
important in the' present problem, which have not been considered 
previously in the theory of curved shells: (l) The effect of transverse 

shear deformation and (2) the effect of transverse normal stress 
deformation. The first of these two effects has been known to be of 
importance in the theory of plates and beams. The second effect was 
found to occur in a manner which is typical for shells and has no 
counterpart in fla 1 >-plate theory. 

The general results of this report have been applied to the 
solution of problems concerning flat plates, circular rings, circular 
cylindrical shells, and spherical shells. In each case numerical 
examples have been given, illustrating the magnitude of the effects 
of transverse shear and normal stress deformation. 

The results of this investigation indicate the necessity of 
taking account of transverse shear and normal stress in sandwich— type 
shells, as soon as there is an order— of-^nagnitude difference between 
the elastic constants of the core layer and of the face layers of the 
composite shell. It was found, that the changes due to transverse 
shear and normal stress deformation in the core may be so large as to 
be no mere corrections to the results of the theory without transverse 
core flexibility. 

The actual magnitude of the changes is greatly dependent on the 
geometry and loading condition of the structure under consideration 
so that no general rules may be given which indicate for which, elastic 
modulus ratio the changes begin to be significant. 

c 

Solutions of problems in the present theory may in general be 
obtained by mathematical methods which are similar to those employed 
in the theory of plates and shells without the effect of transverse 
shear and normal stress deformation included. The present work does 
not include consideration of buckling and finite deflection effects. 



2 


NACA TN No. 1832 


INTRODUCTION 


In this report an extension of the classical theory of small 
"bending and stretching of thin elastic shells is considered. Instead 
of a homogeneous shell, investigation is made of a shell constructed 
in three layers: A core layer of thickness h with elastic 

constants E c , G c , and v c and two face layers of thickness t with 
elastic constants Ef, Gf, and Vf. In the developments certain 

restrictive assumptions are made which somewhat limit the general 
applicability of the results. In so doing formulas are obtained which 
are as compact as possible while still describing the essential 
characteristics of the sandwich-type shell. 

The thickness ratio t/h is assumed Bmall compared with unity; 
at the s ame time the ratio Eft/E c h is assumed large compared with 
unity. This latter assumption means that the face material is so 
much stiffer than the core material that the contribution of the core 
layer to stress couples and tangential stress resultants of the composite 
shell is negligible. It is known that for flat plates these assumptions 
necessitate the taking into account of the effect of transverse shear 
deformation. (See, for instance, reference 1.) Th9 same would be 
expected to be true for curved shells, and the present report, therefore, 
gives a system of equations in which this effect is incorporated. 

A further effect which, it appears, has not been considered 
previously in th9 analysis of small deflections of sandwich structures 
is the effect of transverse normal stress deformation. In the present 
report it is shown that this effect arises in a manner which is typical 
for shells and has no counterpart in plate theory. It may be likened, 
roughly, to what happens in the bending of curved tubes. 

The process by which the general results of this report are 
obtained is as follows: First, each of the face layers of thickness t 

is assumed to behave like a thin shell without bending stiffness. The 
loads applied to these face shells, henceforth called face membranes, 
are- of two kinds: (l) External loads and (2) loads caused by the stresses 

in the core layer. Next, the core layer of thickness h is assumed 
to behave like a three-dimensional elastic continuum in which those 
stresses which are -parallel to the faces are negligible compared with 
the transverse shear and normal stresses. On the basis, of these two 
assumptions three steps are carried out. First, the equilibrium 
equations of the core layer and of the face layers are obtained. Then 
an appropriate expression for the strain energy of the composite 
structure is derived. Finally, Castigliano's theorem of minimum 
complementary energy is used to obtain the relations which connect 
stress resultants and couples of the composite shell with the quantities 
which describe the state of deformation of the composite shell. 
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The system of equations which is obtained in the foregoing manner is 
specialized for the following cases: 

(1) Flat plate 

(2) Circular ring 

(3) Circular cylindrical shell 

(h) Spherical shell with axisymmetrical deformation 

In each case a number of problems are solved explicitly and the 
appreciable effect of transverse shear and/or normal stress deformation 
is illustrated numerically. 

This work was conducted at the Massachusetts Institute of 
Technology under the sponsorship and with the financial assistance 
of the National Advisory Committee for Aeronautics. 

SYMBOLS 


h core— layer thickness 

t face— layer thickness 

curvilinear coordinates on middle surface of 
composite shell 

£ distance coordinate measured along normal to 

middle surface of shell 

a, , a- coefficients of linear element on middle surface 

X of shell 

R. , principal radii of curvature of middle surface of 

shell 

NnTnn direct stress resultants in upper face membrane; 

n = 1,2; m = 1,2 

N riTn | direct stress resultants in lower face membrane 

Pnu-> P nl tangential components of external load intensity 

on upper and lower membranes 

q , q-j normal components of external load intensity on 

upper and lover membranes 
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T i£ ' T 2i 

components of transverse shear stress in core layer 


component of transverse normal stress in core layer 

T n£u J T n£ l 

values of transverse shear stresses for £ = ± h/2; 
n = 1,2 ' 

<V a t,i 

values of transverse normal stresses for £ = — h/2 

T n£m 

values of transverse shear stresses at middle surface 
of shell 

Q 2 

transverse shear* stress resultants 

V 

direct stress resultants, parallel to middle surface 
for composite shell; n = 1,2; m = 1,2 

Mnm 

stress couples for composite shell; n = 1,2; m = 1,2 

Pn 

tangential components of external load intensity for 
composite shell; n = 1,2 

q 

normal component of external load intensity for 
composite shell 

s 

external load intensity term defined hy equation (22) 

Jt 

strain energy 

E f.> G f > v 

elastic moduli of isotropic face-layer material; 
V = V f 

E c G c 

elastic moduli in transverse direction of core— layer 
material 

u li u 2 

effective tangential components of displacement of 
elements of composite shell 

w 

effective normal component of displacement of elements 
of composite shell 

Pi, P2 

effective components of change of slope of normal 
to middle surface of. composite shell 

e £m 

component of strain 

C* = 2tE f 




NACA TN No. 1832 
D* = (1/2) t(h + t)2E f 


C = C*/(l - V 2 ) 
D 

y 

r, 0 
a 

i, 0 

*•1* ^2j ^12 

H = jr/Z 
Z 

Itt^ j Dl2 

k 

0 , 0 
(D 

K. 


bending stiffness factor (D = D*/(l — v 2 )) 

Cartesian coordinates in plane of flat plate 

polar coordinates in plane of flat plate 

radius of circular ring, cylindrical shell, and 
spherical shell 

surface coordinates on cylindrical shell 
parameters defined by equation (63) 

half wave length of sinusoidal load distribution 
quantities defined by equation (197) 

complex quantity defined by equation (200) 
surface coordinates on spherical shell 
quantity defined by equation (74) 
parameter defined by equation (190) 


I — GENERAL THEORY 
Statics of Sandwich-Type Shell 


In order to derive a complete system of equations for the a y n 
composed of face layers and core layers it is necessary first to 
consider separately the statics of the face layers and that of the 
core layer of the shell. Combination of the results obtained for 
these two components of the composite structure must and will lead 
to those differential equations of equilibrium which hold for elements 
of a shell, whether .this shell is of homogeneous or nonhomogeneous 
construction. In addition, however, relations are obtained which 
are characteristic of the sandwich-type shell. 

Coordinate system on shell .- A curvilinear coordinate system 

£) 18 introduced as follows: Let and | 2 fc e coordinates 
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on the middle surface of the composite shell and let £ he the distance 
of a point of the shell from its middle surface, measured along the 
normal to the middle surface. In order that this system of coordinates 
he an orthogonal system, choose the |q, |g curves as lines of curvature 

on the middle surface (in the case of shells of revolution the lines of 
curvature are identical with the meridians and parallels on the middle 
surface) . 

The linear element in the foregoing system of coordinates is of 
the form 



where aq and 09 are the coefficients of the linear element on the 
middle surface and Bq and Eg are the principal radii of curvature 
of the middle surface (see fig. l) . Formulas for the calculation of 
the quantities Oq and B n are contained in texts on differential 
geometry. They are collected, together with other results, in 
reference 2 , which deals with the theory of homogeneous thin shells. 

Statics of face layers .— The face layers are treated as thin shells 
of thickness t and it is assumed that the tending stiffness of these 
thin shells about their own middle surface may be neglected.! Because of 
this neglect from now on they will be designated as face membranes. 

The middle surfaces of the face membranes evidently are given with 
reference to the three-dimensional system of curvilinear coordinates 

by £ = ?j(h + t) and | = — ^(h + t). From equation (l) it follows 

that the linear element on the middle surfaces of the face membranes 
is given by 


ds^ = Ctq' 




dli 2 ;+' a^ 2 


( l 1 nsy) + d ? ! 


(s) 


The components of external load intensity on the upper and lower 
membranes are designated by p lu , p^, and q^. and by , p 1? , p 2 ^, and q ? , 

respectively (fig. 2 ). The core— layer stresses which act on the upper 
and lower membranes are given as Tq£ u , t 2 £ u , and 0 £ u and by ^ r 2 tl> 


^This, of course, means that no local buckling phenomena are 
considered in the present work. 
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and respectively. Finally, the direct stress resultants in the 

upper and lower face membranes are designated by N llu , % 2 u> N 21u , 
and N 22u and by N-^j, N 21Z , and N 222 , respectively (fig. 2). 

There are then three equations of force equilibrium for the 
elements of each of the two membranes. Writing 

/, h + t\ 

°nu " “n^ 1 + 2 E n ) 

- ^f) 

2 

the equations for the upper-face membrane are the following : 



22 gjBU» + - n 22u - ^ 5 ,,) = 0 (4) 


+ N 21u - N llu + a lu a2u(p2u - - 0 (5) 

0^2 O&i o5 2 


a lu a 2u 


H llu + N 22 u 

e i( 1 + 5 st) e 4 + T5t) 


— q u + 


= 0 


(6) 


The corresponding equations for the lowei^-face membrane are 


+ N 12 ; - N 22 1 + a ll“2l (^1! + T l!U) ' 0 ^ 


^2*2i ♦ + N 21I ||p - H m 1^1 + a ll a 2l (P 2I ♦ T 2C ,) - 0 (8) 




2 These are obtained from the corresponding equations of reference 2 
with ctn changed to and with stress couples and transverse shear 

stress resultants omitted. To make up for this omission, the loads on 
the two membranes are assumed to act at their middle surfaces; this means 
terms of the order t/R are neglected (but not terms of order h/R) . 
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/On 




Nil! 


N, 


22 1 


.^-w) 


21 . 


<ll ~ Hi 


= 0 


(9) 


As tending moments and transverse shears are assumed not to he 
acting in the individual membranes the moment equilibrium equations 
become the symmetry relations 


N 12u = W 21u 
W 122 = N 212 


( 10 ) 


Before analyzing the state of stress in the core layer it is 
convenient to see what relations follow from equations (4) to ( 9 ) for 
the composite shell. ■ 

Statics of composite shell .— It may he seen that, in view of the 
fact that all face— parallel stresses in the core layer are neglected, 
the following expressions for the face— parallel stress resultants and 
couples of the composite shell are obtained: 


Nil - (l 

“l2 - (l 
“21 - (l 
“22 - (l 


h + 

2Rr 

£ 

r> 'nu 

+ (> 

h 2Br 
£ 

f>H! 

(11) 

h + 

SR, 

£ 

^)%2U 

4 

h + IAtvt 
^rTJ N 121 

(12) 

h ± 

2E] 



h + 
2B : 


(13) 

h + 
2R] 

^) H 22u 

b 1 

h + 

2 h 

~^) W 222 

(14) 

( 1 + 

2R 2 ) N 11u 

-b 

1 ro «• 

rH 

43 OJ 
> 

(15) 

b 

2B 2 ) N 12u 

-b 

h + tA-w 
2R 2 fl2 1 

( 16 ) 
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«ai -H* 

i 1 + W) Ni2 “ i 1 


( 1 + 2R ± ) N 22u 2B 1 ) N 22iJ 


( 17 ) 

(18) 


In the same way the following expressions are obtained for components 
of external force and moment intensity: 


% 


Pn - ( 

* + ^X 1 + ^sr) pnu + (: 

L W)( 

1 h + tA 

1 -2^77 Pn 2 

( 19 ) 

<1 = 

i 1 + WX 1 + WH + ( : 


(20) 

h + t 
2 

[( 1 + 2E 2 )( 1 + 2 R X ) Pnu I 

L ' ' ' 

4^ ^ 

CVJ 

A 

1 

.H 

C 2®! X* 1 ] 

(21) 


Further , a load term of the following form will be encountered: 


8 - i 


} ♦ * *>£*>. - (* - ^)(l - 


( 22 ) 


which hears a relation to equation (20) similar to that which equation (21) 
hears to equation (19) . This last term would represent, for a homogeneous 
shell, the average transverse normal stress at any station of the shell, 
assum i ng that the loads q u and qj alone are responsible for this stress. 
For a homogeneous isotropic shell this term is of no importance. For a 
sandwich-type shell, as will he seen, it may sometimes he of importance. 

In order to obtain force and moment equilibrium equations for the 
composite shell the face-membrane equilibrium equations (4) to (9) are 
combined suitably. Adding equations ( 4 ) and (7), and (5) and (8), 
respectively, the two equilibrium equations for the force components 
parallel to the middle surface of the shell are obtained. In order to 
reduce them to known form (see reference 2) the following relations are 

used between the core-layer-surface shear stresses r c and t » 

nQu n<->2 

and the transverse shear stress resultants Q-^ and Qg. 
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h + t 

2Ei 



h + t\ 

2e 2 ; 


T 


n£u 



3ft 

*n 


(23) 


h + t 


TH^ 1 + TiB^Vu + (l - ^(l - 


« On (24)* 


Equationa (23) and (24) will subsequently be shown to be in agreement 
with the usual definition for the transverse shear stress resultants by 
consideration of the stress distribution of the core layer. 

With equations (23) and (24) , there are obtained by combination of 
equations (4) and (7), and (5) and (8) - carrying out addition as well 
as subtraction — the following four equations: 


2S&1 + 32J&L + n 12 |l - + 


86 - 


8|/ 


+ n, 


= 0 


SotgN^g 8(X^N 22 

86-1 860 


+ N, 


21 


Sag 

817 


- N 


11 


8a^ 

W c 


9 


a-^cxp 


P 2 


= 0 


r ■ 

8g^i + 8a^i + % 2 ~ **22 ZF + ~ %) = 0 


dig 


8 a 2 Mi2 8aq^ 2 


8a 2 


8 aq 


86 -, 


86 , 


+ ~ M 11 XI7 + a l a 2( m 2 " %>) = 0 


86 


(25) * 

(26) * 

(27) * 

(28) * 


Two further equations are obtained by adding and subtracting, 
respectively, equations (6) and (9). Adding equations (6) and (9) and 
taking account of equations (ll), (l4), and (20), there follows: 



( 29 ) 
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In order that this reduces to the correct equation of transverse force 
equilibrium as given in reference 2, 



Equation (30), just as equations (23) and ( 2 k) , can again be verified 
independently by consideration of the state of stress in the core layer. 
On the basis of equation (30), equation (29) is written in the form . 


^FT 



+ a^agq = 0 


(3D* 


The last equation, use of which is required for the sandwich-type 
shell and which has not previously been given, is obtained by 
subtracting equation (9) from equation (6). Taking account of 
equations (15), (l8), and (22), there results 


+ (1 - wX 1 - wH - 0 


1 + w) o?u 


(32) 


Provisionally, there is written 


C 1 * WX 1 + ^>6. * i 1 - ^fX 1 ~ ‘ (33) 
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and it will subsequently b9 shown that CT£ m represents the value 
of <J£ at the middle surface of the shell. Combining equations (33) 
and (32) yields 


+ 


1 /%! 
h + typi 



- s = 


0 


( 3«0 


Equation ( 34 ) has no relation to the sixth equation of equilibrium 
for an element of the shell which expresses the condition of moment 
equilibrium about the normal to the middle surface. That equation which, 
as is known, is an identity when resultants and couples are expressed 
in terms of stresses does not occur in the present derivations, or rather 
it is contained in equations (12), (13), (l 6), and (17), which give 
explicitly the slight differences between N^g ' and. and Mto 

and . 

Stress distribution in core layer .— In order to verify independently 
equations (23), ( 2 k), and (30), as well as for the subsequent derivation 
of appropriate stress-strain relations, it is necessary to determine the 
distribution of stress in the core layer. 

Assuming that the components of stress a^, Cg, and t-^ i n the 

core which would contribute to stress resultants and couples of the 
composite shell are of negligible importance, 3 these components of 
stress may be set equal to zero and only the components of transverse 
shear stress and transverse normal stress t^, tv >£, and cr£ may be 

retained. The differential equations of equilibrium for these three 
remaining components of stress in the system of curvilinear coordinates 
defined by equation (l) are obtained, from the general form of these 
differential equations in reference 3, in the following form: 



3 lt is for this purpose that the order— of-magnitude relation 
hEp^tEj. « 1 is assumed. 
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The values of the three stress components at the middle 
surface ( £ = 0) are designated by the subscript m. Integration of 
equations (35) to (37) then gives 


r lt 


T lU 

(1 + C/Ri) 2 (i + tfa>) 


(38) 



T g£m 

(1 + t/R 2 ) 2 (i + £/%) 



(39) 


m 


The transverse shear stress resultants Qj 
from equations (38) and (39) in the form. 1 * 


and Q2 are obtained 


On = 


fUh+t)/2 


u 


l-(h+t)/2 



^n T n£m/ r 

V 




(h + t)T^ 



( 41 ) 


^The integration must be extended over the thickness of the core 
layer and also over half the thickness of the face layers, in accordance 
with the prior assumption that the stresses T n £ u , cr£ u , 

and 0^1 may be taken to act at the middle surfaces of the respective 
face membranes. 
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Now, as intended, the proof is carried out of equations (23), (24), 
and (30), which were used to obtain the differential equations for the 
composite shell. 

To verify equation (23), from equations (38) and (39) for the 
left-hand side of equation (23), the following equation is obtained: 


T nU 


T n£m 


1 I h ± t _ h + t 


2^ 


2*n 


h + t T nlm 



and this, in conjunction with equation (4l), verifies equation (23). 

To verify equation (24) in the same manner, from equations (38) 
and (39) for the left side of equation (24), the following equation 
is obtained: 


h + t / T nCm 


h + t 
2E n 



0> + t) T n£m 



and this, in conjunction with equation (4l), verifies equation (24). 

To verify equation (30), equation (40) is used to write for the 
left side of equation (30) 


d /°2 t lSm \ + S I a l T 2£m ^ + d (jVjU \ + d / a l T 2U 


dijL . h + t I dlp| n . h + t ) (, h + t I dlp|, h + t 


2Ri 


■2i 


2Kr 


• os l 1 - 


2R-, 


’2 1 - 


2Rr 


= d / 2a g T lU \ , d / 2a l T 2tm 


^lL _ /h + t\ 2 1 ^ 2| 

2R-, 


1 - 


h + t 

2Rr 


and this, in conjunction with equation (4l), verifies equation (30) 
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The section on the stress distribution in the core layer is 
concluded by listing the form which equations ( 38 ) to (40) for the 
stresses in the core layers assume for "thin" shells, that is, for 
shells for which h/R « 1. From equations ( 38 ) and (39), in 
conjunction with equation (4l), it follows that 


Q 1 

T i<; - ftt 

Q 2 

T 25 = h-Pt 


(hs)* 


From equation (40), in conjunction with equation (4l), it follows that 


(XiCCrjCF £ 


“l^U - 


h + t^a Sl 



(W 


It is necessary to note for some of the following considerations 
that, in view of equation ( 31 ) , instead of equation (^ 3 a) there may 
be written - 


♦§*-«) (“3,). 


It is seen that in this approximation the transverse shear stresses are 
uniform across the thickness of the core layer, while the transverse 
normal stress is composed of two terms, one uniform across the thickness 
and the other varying linearly across the thickness. 

Wo further calculations are needed with reference to the state 
of stress in the composite shell. The next step is to complete the 
system of differential equations for stress resultants and couples by 
deriving an appropriate system of stress— strain relations. 


Strain Energy of Sandwich— Type Shell 

In calculating the strain energy of face membranes and core layer 
it is assumed that both are isotropic and elastic, with elastic 
constants E f , v f = V, G f = Ef/2(l + v) and E c , V c , G c = E c /2(l + V c ). 



1 6 


NACA TN Wo. 1832 


Poisson's ratio for the face membranes is written without a subscript , 
because, in view of the assumed stress distribution, there is no explicit 
occurrence of Poisson's ratio v c for the core layer. 

The strain energy for the composite shell is the sum of the strain 
energies for the face membranes and for the core layer 


it = Jtf + ir c (44) 

For the purpose of obtaining stress-strain relations, both jtf and n c 
are expressed in terms of stresses rather than in terms of strains. 

Strain energy of face layers .— Considering that the element of 
area on the middle surfaces of the membranes is of the 

form ± ^ ^1 ± h gR <^1 dl 2 and that the stresses 

in the membranes are the stress resultants divided by the membrane 
thickness t, there is, from well-known principles, the following 
relation: 

”f - | jf ^["llu 2 + ”22u 2 - + 2(1 + v)N 12u 2 ] 

* (! + + «1 

+ 2 Jf t%[ B lll 2 + N 22l 2 - 2VB 11I B 221 + 2 <1 * v)N 12 , 2 ] 

* ( r - - ^) a l a 2 «1 «2 

Equation (45) is transformed into an expression containing stress 
resultants and couples of the composite shell by mean a of equations (ll) 
to (l8) which lead to the relations 


1 
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2 ( 1 + ^~) N Hu = N u + FTt 


2 ( 1 “ = N n - rrr 

with corresponding formulas for N-|p and N 22 .5 

In this transformation the cases are limited to those for 
which h/R « 1. Then, with the two constants C* and D* defined by 


**L1 


%1 


(46) 


C* = 2tE f 

D* = 1 t(h + t)?E f 




(47) 


the following expression for is obtained: 


«t =2 


uu 


c* 


[nh 2 + Npg 2 - 2vN n N 22 + 2(1 + v)N 12 2 ] 


^[Wll 2 + “22 2 - 2v%l>fe- + 2(1 + V)^]} a-jOjj dl x dft 2 (48)- 


It may he remarked that equation (48) could have been given directly, by 
analogy with known results for the isotropic homogeneous shell. 

Strain energy of core layer .— With the stresses a-^, a 2 > and t^ 2 
assumed to vanish, there results for the strain energy of the core layer 


^Note that equations (46) and corresponding equations can be used 
to calculate the stresses in the two different face membranes, once 
stress resultants and couples in the composite shell are known. 



18 


NACA TN No. 1832 



Again the terms £ /R compared with unity are neglected and, consistent 
with this neglect, the values of the stresses and cr£ are taken 

from equations (42) and (43) . 

The value of cr£ may he chosen from either, equation (43a) or 
equation (43b) . The form of the results depends somewhat on which of 
the two equations is chosen, in the sense that the meaning of the 
deformation quantities which are to he determined depends on which 
of the two equations is taken. This question is decided in the 
following manner: As all resultants and couples enter the expression 

for the strain energy only as themselves and not in differentiated 
form, except when equation (43a) is used, the selection of equation (43b) 
for C£ is proposed, thereby excluding derivatives of stress resultants 

and couples from the expression for the strain energy jt. 


Introducing then equation (43b) into equation (49) yields 
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(5.0) 


The integration with respect to £ is carried out and 
equation (50) becomes 


r\ 
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Q l 2 + ^ 2 , h + t 
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(h + t)G c + E c 
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L J J 


o^ag d^ dl 2 (51) 


It was to be expected that the terms containing the modulus of 
rigidity G c would occur in the foregoing form. The contribution of 
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the present report up to this point, besides giving the new equation (34) 
for cr£ m , is thought to be the determination of the form in which the 

effect of transverse normal stress deformability manifests itself in 
the strain energy of the sandwich shell. 


Stress— Strain Relations for Composite Shell 

In what follows a system of stress— strain relations for the com- 
posite shell is obtained by the use of Castigliano’s theorem of minimum 
complementary energy. The manner in which the theorem is used here 
appears to have been employed first by E. Trefftz (reference 4) for the 
purpose of avoiding geometrical considerations in the derivation of the 
stress— strain relations for thin homogeneous shells with small defor- 
mations, without consideration of the effects of transverse shear and 
normal stress deformation. 

Assuming for the present purpose that all boundary conditions for 
the shell under consideration are stress conditions, the theorem 
consists in the statement that among all statically correct states of 
stress the actually occurring state of stress makes the strain energy 
of the system a minimum. In the application of the theorem the fact 
is taken into account that statically correct states of stress only 
are to be compared, by means of the Lagrangian multiplier method. 

Before minimizing it an integral is added to it which contains the 
six equilibrium equations (25) to (28), ( 31 ), and (34),' each of the 
six equations multiplied by a lagrangian multiplier. It can then 
be shown, by using Castigliano’s theorem with prescribed boundary 
displacements instead of with prescribed boundary stresses, that 
each of the six multipliers has the meaning of one of the displacement 
quantities which occur in the shell problem . 6 

With the foregoing understanding of the meaning of the multipliers, 
the multiplier of equation (25) is designated by ujj that of 

equation (2 6), by ugj that of equation (27), by P^j that of 

equation (28), by P 2 ; that of equation (31), by w; and finally that 

of equation (34), by k. It is known that up, u 2 , and w represent 
the effective components of displacement in the ij, ^ 2 , and £ direc- 
tions, respectively. Further, it is known that Pi and p 2 represent 
the angles through which the normal ±0 the middle surface of the shell 
turns toward the £1 and I 2 curves, respectively. There is no 

^For the special case of the flat plate this has been carried out 
explicitly in reference 1. For the case of the homogeneous shell, 
without effect of transverse shear and normal stress deformation, the 
proof has been given in reference 4. The proof for the more general 
case which is here considered is not included as it does not offer any 
clearer insight into the problem and tends to lengthen the analytical 
discussion. 
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Immediate simple geometrical interpretation for k and, while such 
interpretation in terms of an average transverse normal strain might 
be deduced herein, k is considered as an auxiliary variable presently 
to be eliminated. 

Combining now equations (*<4), (48), (51), and (25) to (28), (31), 
and (34) in the manner indicated, the following variational equation 
results: 
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(52) 
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The variations in equation ( 52 ) are carried out and integration is 
done "by parts to eliminate derivatives of variations in the double 
integral. The line integrals along the boundary which occur due to 
this integration by parts vanish, because it has been assumed that, all 
stresses are prescribed at the boundary and therefore their variations 
vanish at the boundary. 


The resultant variational equation is 
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( 53 ) 
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As all nine variations in equation (53) are independent of each other, 
it follows that the contents of all nine "brackets in equation (53) must 
vanish separately. Thus the following nine stress— strain relations are 
obtained for the sandwich shell, indicating with an asterisk those 
which appear in final form. 



(54) 


N22/ (h + t)cA _ Ny / _ (h + t)C»N 
C* V 12E c E 2 2 J C* \ 12 E cEl R 2 J 

_ jl_ Me. U 1 dag w_ (h + t)q 
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D * «1 Sli °i a 2 dig “l 0 ^ (h + t)E x 

1^2 - vM ii _ _i_ Me + ^1 Me _ k 1 

D* 02 dig c^ag ail c^ag (h + t)E 2 


(57) 


( 58 ) 
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% 

(h + t)G c 


= 3, + -L ik_ _ lil 
P1 a l Mi El 


% •_ Bo . _1_ dw ^2 

(h + t)G c P2 02 5I£ E 2 


jk = k_ . 1 . 

E c ai«i2 (h + t) 


(60)* 

(61)* 
(62) , 


It may "be verified that the meaning of the quantities u 2 , w, fii, 

and 0 2 is as haa "been indicated by comparing equations (54) to (6l) 

with the corresponding equations of reference 2 for the homogeneous 
shell with E c = G c = 00. 

The system of equations (54) to (62) may be brought into a 
slightly more concise form as follows: Define the quantities X-^ X 2 , 

and X-^ 2 by 


. 1 (h + t)t E f 

1 " 2 El 2 E c 



(h + t)t E^ 
~R? E c 




(63) 


X12 


1 (h + t)t 

2 E lE2 E c 


and eliminate k from equations (57) and (58) by means of equation (62) 
and the equilibrium equation (34). Retain equations (5 6) and (59) to (6l) 
in the foregoing form and write for equations (54) and (55) 


(1 * f *l)«U - (v - \ te)fe> . c(h- |i + 


u 2. Sai v\ 

Sly 


(h + t)C* 
12E cEl 


( 64 )* 
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(* ♦ H - (* - * - c *(4 n % ♦ *) 

+ (h + t)C* 

Equations ( 57 ) and (58) 150001110, 
according to equation ( 47 ), 


12 E c R 2 

if D*/(h + t)E c = | t(h + t)E f /E c , 


(1 - Xi)Mu - (v - Xi2 )^ 2 = D.^. ||1 . | 2 i^ + s ( 66 )* 

(1 + X2)Mj 2 - ( v - X 12) M 11 = ||| + + 8 (67) * 


With these, last transformations there is obtained a system of 
equations which is formally equivalent to the corresponding system of 
equations for the homogeneous shell. The 5 equilibrium equations ( 25 ). 
to (28) and (31) and the 8 stress— strain relations (56), (59), (60), (6l), 
and ( 64 ) to (67) are used for the determination of 13 quantities: Five 

stress resultants N^, N 22 , N^, Qj_, and Qq; three stress couples M-q, 

M 22 , and and five displacements and changes of slope u^, Ug, w, p ^ , 

and P 2 . The quantity cx^ which occurs in. the sixth equilibrium 
equation (equation ( 34 )) may be determined directly, once the shell 
bending and stretching problem has been solved. 

It is seen that the effect of transverse sheen* deformation enters 
equations ( 60 ) and (6l) only and that, when G c = <*>, these equations 
give the values of the known theory of homogeneous shells without 
transverse shear deformation (references 2, 3 , and 4 ). 

The effect of transverse normal stress deformation enters equa- 
tions ( 64 ) to (67) only. It is seen that it is, in part, responsible 
for the occurrence of apparent stiffness factors C*/(l + X) 
and D*/(l + X) . Thus, according to equation (63), the effect of 
finite E c is to make the shell more flexible in bending and stretching 
them it would be with E c = », This effect, however, is present only 
in curved structures and not in plates and straight beams, as the 
quantities X have one or both of the radii of curvature in the 
denominator. A further effect of finite E c is occurrence of the 
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external load terms q and s in the stress— strain relations. Both 
these effects represent , roughly speaking, what happens to the shape 
of an element of the composite shell if the length of the core fibers 
in transverse direction is changed, without any stretching or 
compressing of the face-membrane elements. 

Having derived the general system of equations for the small 
bending and stretching of sandwich-type shells, it remains to apply 
these equations to specific problems which may be of interest and to 
determine the quantitative effect of the terms which are characteristic 
of the sandwich-type shell. Some of this work is done in part II of 
the present report, which follows. 

It may be stated once more that for these specific applications 
the five equilibrium equations (25) to (28) and (31) and the eight 
stress-strain relations (56), (59), (60), (6l), and (64) to (67) are 
used. 


II - APPLICATIONS OF GENERAL THEORY 
Flat Plates 


The problem of the flat plate is considered first in order to show 
that the results of reference 1 are contained in the present results 
and in order to solve some problems in the theory of plates which have 
not been solved in reference 1. 

Rectangular plates .— Using notation which is customary in plate 
theory there is set 
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The equilibrium equations (25 ) to (28) and (31) become 


8N t 3N- 


+ + Pz ° 

SNxy 6N V 

if + + 


(69) 


§f + ;v +,I = 0 


SMrv SMn- 

^ + 5 ?-^ + “V’° 


(70) 


The stress— strain relations (56), (59)* (60), (6l), and (64) 
to (67) become 


N x - VNy = C* 


8x 


»y - *«z - 0* ^ 


2(1 ♦ v)^ . C(| ♦ £) 


(71) 
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Mj — VMp. = D* 


- vMj = D* 


2(1 + v)M*y = 


As In the small-deflection theory of homogeneous plates, the 
equations for stretching (equations (69) and (71)) are independent of 
the remaining equations for transverse bending. Equations (69) and (71) 
for the stretching are not affected by the elastic properties of the 
core layers . 

Equations (JO) , (72), and (73) have been treated in reference 1 by 
means of a stress function ty, which, together with the deflection w, 
was taken as one of two basic variables. In what follows an alternate 
treatment is given, in which the problem is reduced to three simultaneous 
equations for the quantities 3 X , 3 y , and w. On the basis of these three 
simultaneous equations a problem not considered in reference 1 is treated, 
namely, the bending of a rectangular plate which is simply supported on 
all four edges. 7 

To reduce equations (70), (72), and (73) to three simultaneous 
equations for P x , 3y, and w, first a quantity 0) is defined by 



= 83. 


dy 


( 74 ) 


Introducing equation (72) into the first of equations (70), in view of 
equation (74), there is obtained 


to + V2w = - q/(h + t)G c (75) 

\ 

Next, Qj, Mj-, and are taken from equations (72) and (73) and the 

result is substituted in the second of equations (70). This gives, after 
slight transformations. 


7 

This same problem has also been solved by L. H. Donnell by a 
method which differs from the one employed here. (See reference 5 
where the case of the homogeneous plate is considered.) 
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rhj v2p * - 2(h + t)G c p x + ^[rrr - 2 < h + t)G c v ] + % = o (76) 

In an analogous manner the following further equation is obtained: 

jB-Tv - 2(h + t,G cPy + - 2 <» + tJv] + my = 0 (77) 


In order to solve equations ( 75 ) to ( 77 ) two equations are next obtained 
involving w and to only. Differentiating equation (76) with respect 
to x and equation (77) with respect to y and adding the two resultant 
equations , in view of equation (7^0, gives 


^^2 V2 “ " 2(h + t)G c( <D + ^ w ) + ^ » 0 


and, making use of equation (75 ) , 



The following procedure may now be carried out: (a) Solve 

equation (78) for u>, (b) with this value of <0 solve equation (75) 
for w, (c) substitute oi and w in equations (76) and (77) and solve 
for 3 X and 3y, and (d) eliminate extraneous terms in 3 X and 3 y 
by considering equation (7^) • 

Before deriving the solution of a problem along these lines, 
the explicit differential equation for w which follows by combining 
equations (75) and (78) may be given 


V^7 2 w 



_ . A 


(h + t)G c 


( 79 ) 


Note that the effect of transverse shear occurs on the right side of 
the equation only. In order to compare the magnitude of the q terms 
on .the right of equation ( 79 ), assume that relevant changes of q 
occur over distances of order Z (where Z may or may not be a 
representative diameter of the plate). Then, as ordeiwjf— magnitude 
relations, there results 
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From equation (80), it follows that transverse shear ceases to he a 
secondary effect as soon as 2 is of order \/ht or of smaller 

order . 

Bending of rectangular plate with simply supported edges .— The edges 
of the plate are assumed to he at x = 0 ,a and y = 0,h and along these 
edges moments and deflections are assumed to vanish. Further, 


m x = my = 0 


* = X 8in Bin ^y 

m=l n=l 


(81) 


where 


= m «/a 

Mn = nn/b 


(82) 


From equation (78), it follows that 





, 2 2 
V + -Mn 


sin sin uny + 0% 


(83) 


where is a harmonic function. Putting equation (83) into 

equation (75), 
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Imn 


1/D 
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Mn 2 (h + t)G c 


sin \jx. sin i^y - % 


which is integrated to 


"ill 


± mn 


(v 2 


1 + 


DU. 


Mn 


m + jjn_ 
(h + t)G c 


sin XjjjX sin ^y + w h (84) 


where w h is the general solution of V 2 w } 1 = - 0 %. It is to he expected 
and may he shown explicitly that for the plate which is simply supported 
all around w h = = 0 and, as in the Navier solution for the plate 

without transverse shear deformation, the particular integral is the 
complete solution of the problem. 

Equation (84) may he rewritten in the more explicit form 


w 



< lmn'j 

[ ; £ _ 

Ef (h + t)l 

-[m 2 + n 2 (a 2 /h 2 )J ' 

L 2 (1 ■ 

- V 2 )G c a 2 


m 2 + n 2 (a 2 /h 2 )J 
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X sin 2^ x sin ^ y 
a h J 


(85) 


When G c = «>, equation ( 85 ) reduces to Navier’ s solution Wfj. Equation ( 85 ) 
is more readily interpreted hy means of the ratio w/w K of deflection 

with and without transverse shear deformation. On the basis of 
equation ( 85 ), there may he obtained the following equation ( 86 ),® 


^Setting *2/2 (l - V 2 ) = 5.4 and (E f /G c )t(h + t)/a 2 = p, 

equation (86) takes on a form which contains as a special case the 
result of equation (l8) of reference 5 . 
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( 86 ) 


For the case of a uniform load intensity q = Constant and for the 
center of the plate (x = a/2, y = h/2) equation (86) "becomes 


JL = 2. + — ^f (h + t)t 

w W ' 2 (! - v 2 )G c a 2 



sin mjr/2 sin ng/2 


mn 


m 2 + n 2 (a 2 /h 2 )] 
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[m^ + n 2 (a 2 /b 2 )J 


(87) 


The ratio of the series is I.98, when a/b = 1, and the ratio of the 
series is 1.11 when a/b = l/2. 

For the case of a concentrated load at the center of the plate 
the deflection ratio at the point of load application assumes the form 


(sin mrt/2 sin nm/2) 2 

m 2 + n 2 (a 2 /h 2 ) (g8 

(sin mm/2 sin nm/2) 2 
[m 2 + n^ (a 2 /h2 )J 2 

Wow it is easily shown that the numerator series in equation (88) 
does not converge and consequently w/wn = «> in this case. A more 
detailed consideration shows that in any plate theory which takes 
transverse shear deformation into account the deflection under the point 
of application of a concentrated load must "become infinite in contrast 
to what happens when transverse shear deformation is not taken into 
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account. This difference, of course, vanishes as soon as the load 
intensity becomes finite, and then the theory with transverse shear 
deformation taken into account is more accurate than the theory which 
does not take into account this effect. 

For the sake of numerical illustration take again the square 
plate (a/b = l) with uniform load distribution. According to equation (87) 
the deflection at the center is increased because of transverse shear by 
the factor 


- 1 + 9.7 ft 

wu G c 


(87a) 


Take h = 1.0 inch, t = 0.1 inch, a = 10 inches. Ef/G c = 200, 
and V = l/3. Then, according to equation (87a), w/wr =1 + 2.3, so 
that in this case the deflection with transverse shear is more than three 
times the deflection when shear deformation in the core is neglected. 

Returning now to equation (84) for w and equation (83) for <u 
and substituting these two equations in equations (7 6 ) and (77) in order 
to determine the changes of slope 0 X and 0y, after slight transf ord- 
inations there results 


Px = 





Qmn^m. 

(^m^ + Mn 2 ) 


2 cos X-mX sin p n y 


Bin ^ CO0 ^y 

(*m 2 + hn 2 ) 


(89) 


Equations (89) are remarkable for the reason that they are not affected 
by transverse shear' deformability. According to equations (73)* the 
same is then true of the bending and twisting couples M*, My, and M^y. 

It is not easy to see why, in this statically indeterminate problem, the 
magnitude of the internal forces, as well as that of the deflections, does 
not depend on the elastic properties of the core. The analysis, however, 
shows that the distributions of M^, My, and M^y, and therewith 
of Qx and Qy, remain the same as those obtained under the assumption 
that G c =00. in this connection the following remark may be made. 
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Evidently the following three boundary conditions, v = = p y = 0 

along the edges x = 0,a, have been satisfied. In order that the last 
of these three conditions be satisfied there are necessarily nonvanishing 
edge values of the twisting couples M^. The same is true in the theory 

without transverse shear deformation, where, however, no alternative 
possibility exists, as in that theory only the boundary condi- 
tions w = = 0 are relevant. For the present system of equations 

three boundary conditions must be formulated for every plate edge. Thus 
it is possible although mathematically complicated to solve the problem 
of the rectangular simply supported plate with the edge condition Py = 0 
replaced by the condition Mgy =0. In that case, which will not be 
pursued here, there evidently will be a distribution of internal stresses 
which is modified by the effect of transverse shear deformation. 

Cylindrical bending of plates .- As a further relatively simple 
example of application of equations (70), (72), and (73) problems are 
considered for which 


d( )/dy = 0 

d( )/5x = d( )/dx =( )• 
Mxy = Qy = m.y = 0y = 0 

% = 


and where consequently the problem reduces to the following system of 
equations: 9 

Qj’ + q = 0 

Mg’ — Qrg. + lOj = 0 

Qj = (h + t)G c (p x + w‘) 

(l - V 2 )^ = D*f3 x ' 

^Note that in order to oh tain the problem of the sandwich beam from 
equations ( 9 ^) 8-ud (91) the only changes which are necessary amount to 
setting v = 0 in equations (91). 
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To set into evidence the effect of finite values of G c in 

equation (91 ) , the following system of equations is deduced from 
equation (91): 


Dw v 


(IV) , Dq 

' = q + %’ : — - 


t t 


(h + t)G c 


Mt = - Dw" - 


Dq 


(h + t)G c 




3t - - w* + 




(h + t)G c 


. (92) 


(93) 

(94) 

(95) 


Solutions to the following problems are listed: 

(l) Simply supported plate of span Z carrying a 
load q = qn cos jtx/Z . Boundary conditions: w(±z/2) = U^(±l/ 2) = 0. 


v 




E f 

(h + t)t 

cos nx/Z 

2 (1 

- v 2 ^ 

Z 2 

(n/X) k 


(96) 10 


As the problem is statically determinate as far as moment and force are 
concerned there is no modification of M* and Q z due to the finite 
value of G c . 


10 The factor in brackets may again be written in the form 1 +5.k$ } 

with 0 = (Ef/G c )^[(h + t)t/Z 2 ], using the notation suggested in 
reference 5. 
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(2) Simply supported plate of Brian Z carrying a uniform 
load q = q Q . 


= 3^J 

D16 ^ 


24 


W’ 1 ]- 


f 1 + - 

H (h + t)G c 2 2 . 
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(97) 


From this for the center deflection. 


w(0) = -J- + Q If Xk_ f , t)t ] 

384 D L 3 * 5 (1 - v 2 )g c * 2 J 


(98) 


It is seen that the correction factor for the center deflection is almost 
the same as that for the cosine load curve (equation (96)), the only 
difference being a change of the factor 1^/2 = 4.93 into 24/5 = 4.80 
that is, a reduction of the shear correction factor by at most 3 percent 
is present.il 

(3) Built-in plate of span l carrying a uniform load q = q Q . 

The boundary conditions are: w(±l/2) = 0 Y (±Z/2) = 0 ( and 

not w*(±2/2) =0). 
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1 + — £4d ; 

(h + t)G c 2 2 



(99) 


From this there follows for the center deflection. 


v(0) 


io l 


4 r 


384 D 


1 + 24 


Ef 
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] 


( 100 ) 


Comparison of equations (100) and (98) shows that for the built-in 
plate the effect of transverse shear deformation is very much more 


Note that according to equation (87a) the shear correction factor 
for the square plate of width a = l is more than twice as large as the 
shear correction factor for the plate strip of width 2. 
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pronounced than it is for the simply supported plate, a factor 2 k /5 
in the latter case being replaced hy a factor 2 k in the former case. 12 
As a further result in this problem of the built-in plate, by putting 
equation (99) into equation (93), it is found that the moment 
function does not contain any terms depending bn the effect of 

transverse shear deformation. This again is somewhat surprising as 
in this' case it is not possible to determine the moment function by 

statics alone. -*-3 

Circular plates; rotational symmetry .- As no examples of solutions 
of circular sandwich— plate problems have as yet been published and as 
it is of some interest to determine in which way the shear correction 
factors change in going from a problem for the plate strip to the 
corresponding problem for the circular plate, the equations for 
axisymmetrical transverse bending of circular plates are briefly 
discussed. 

Polar coordinates r ,0 are introduced and notation which is 
customary in plate theory is used. As a consequence of equa- 
tions ( 70 ), (72), and (73), the following system of equations is 
obtained: 


drQj. 

-dF" + r< * 


0 


drMp 

1 — - %-rQr+nnp-O 


( 101 ) 


Qp = (h + t)G c (p r + dw/dr) 


( 102 ) 


Ip 

A somewhat similar percentage increase must take place in going 
from equation ( 86 ) for the rectangular plate with all four edges simply 
supported to a formula (which has not yet been derived) for the rectan- 
gular plate with all four edges built in. 

1 3 as a problem where the moment distribution is in fact dependent 
on the effect of transverse shear there may be mentioned the problem of 
the cylindrically bent plate with both ends built in, which carries a 
load q^ = q^x instead of the load qg = qg. This problem also may be 

solved by means of equations (92) to ( 95 ). 
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Mr — VM 9 = D* d(3 r /dr 

> 

% - vMr = D^Pp/r 


(103) 


According to equation (79), the equation for the deflection w 
will he 




11 drmr _ y 2 q 

q 2 r dr. (h + t)G c 


(104) 


where = (l/r)d[r d( ) /dr] /dr . Having found w by means of 
equation (104), p r may he determined from 


Pr 


= _ &£ + Or 

** (h + t)G c 


. dw 
dr 


+ 



rq dr 


(105) 


and therewith Mr and Me are obtained from equations (IO 3 ). 

In the present problem it seems to he somewhat more convenient to 
proceed as follows: Combine equations (lOl) and (IO 3 ) to obtain as 

equation for the change of slope P r , 


D 1 
r 


d_ 

dr 



= 1 - 


1 dr%- 

r dr 


(106) 


Having Pj., Mr ^ M© are found from equation (IO 3 ) and Qr, from the 
second of equations ( 101 ), 




% - d( 1 3 r + v 


( 107 ) 
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(108) 


Finally, with this value of Qj., w is found "by integrating equation (102), 


w = - / 


Pr dr + 


D 


(h + t)G ( 


1 drPj.\ JW 


(1 tefiA 

\r dr / 


(h + t)G c 


(109) 


Deflection of circular plate with built— in edge .— The tending is 

now considered of a plate with transverse load q = qn(r/a) n and 
with in*. = 0. First, from equation (106), 


D(3 t 


= Cl - + 1 + Co £ loge 21 + p(|) n+3 (110) 

1 a 2 W 3 a °® a ( n + 4)( n + 2 )2\elJ 


Attention is restricted to complete plates with no concentrated 
load at the center, and consequently it is necessary to set c 2 = c^ 

in equation (110) . This gives 


= 0 


D Pr = °1 I + 


Jn2: 


(n + 4)(n + 2) 2Va ' 


n+3 


(110a) 


Putting equation (’110a) into equation (109), there results for the 
transverse deflection w 


D - = - 




< ln a ~ 


© n+ * ♦ H 


(n + 4) (n + 2) 2 


D 


_ i2ct + , 

a 2 (h + t)G c L ( 


In a 3 /r\ n+g l 

(n + 2) 2 '' 8 '^ J 


(111) 


f 
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Taking the case of a plate of radius a with "built— in edge, that 
is, with the "boundary conditions 


^(a) = w(a) = 0 . ( 112 ) 


there results 


D3 r = 


9n a ~ 


(n + 4)(n + 2) 2 




(113) 


and 


Dw = 


_ 

(n + 2)2 


(r/a) n+1 * - 1 _ 1 (r/a ) 2 — 1 

(n + 4)2 2 . n + 4 


a 2 (h + t)G c 


[or 2 - ] 


(114) 


From equation (ll4) there follows for the deflection at the center of 
the plate 


Dw(0) = 


9n a 


2 (n + 2 )(n + 4 ) 2 


1 + 


(n + 4) g 
n + 2 


E f 

(1 - v 2 )g c 


(h + t)t 


(115) 


Consider the following special cases: 

(l) Uniform load distribution q n = qQ. From equation ( 115 ), it 

follows for the ratio of deflection with and without transverse shear 
deformation that 


•w(O) = i + 8 E f (h + t)t . 

[WO)] " (l - v2) Gc a2 

G c =°° 


(116) 
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Equation (ll6) may "be compared with equation (100) for the deflection 
of the infinite plate strip of width 2 with huilt— in edges. 

Setting 2 = 2a, it is seen that, while the transverse shear correction 
factor for the strip has a value 6, the corresponding factor for the 
circular plate is 8. This is consistent with the earlier comparison 
■between the simply supported strip and the simply supported square 
plate, except that there the change is from 4.8 to 9.7. 

(2) Linearly increasing load distribution q = q^r/a. From 
equation (115), it follows that 


v(0) =i + 25 E f (h + t)t 

[>(°)]g c =oo 3 (i-v 2 )g c a2 


(117) 


showing that the correction effect is only slightly greater than in the 
case of the uniform load distribution. 


(3) Load increasing linearly from edge to center, q = q 0 + qq(r/a). 
(ll = ~ goj* From equation (115), it follows by superposition that 


Dw(0) = 


Tpa 

2 X 2 X 16 


<io a 


1 + 8 


E, 


f (h 


(l -v 2 )G c 


+ t)t~ ] 

J 


2 X 3 X 25 


1+22 Ef (h + t)t 

a 2 


3 (l - v 2 )g c ac 


Dw(0) = 


- ^o 1 


,4 r 


32 x 150 


l + 


3000 Sf (h + t)t~ | 
387 (l — V 2 )g c a2 J 


(118) 


Comparing the factor 3OOO/387 = 7-76 which occurs in equation (ll8) 
with the corresponding factors 8 and 8.33 in equations (ll6) and (117) 
it is seen that, in the foregoing three problems at least, there is 
little difference between the transverse shear stress correction factors 
in the case of three different loading conditions for the circular, 
clamped— edge plate. The fact that this agreement should not 'be expected 
to hold generally follows again by considering the case of a point 
load at the center of the plate, for which the shear correction factor 
would again be infinite. 
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The examples of this section should he augmented by the solution 
for the circular plate of radius a, which carries a load distributed 
uniformly over a smaller circle which is concentric with the boundary 
of the plate. 


Circular Bings 

As the simplest example of a curved sandwich structure there are 
considered in this section stresses and deformations of circular rings 
in their own plane. As was found in the general developments of- 
part I of this report, in a curved sandwich structure there will be 
the effect of both transverse shear and normal stress deformation. 1^ 

There are set for the relevant coordinates and variables 


tq = a0 

aq = 1 

Bq = a 

uq = v 

Pi = P 

h* 

H 

II 

II 

<£> 

s. 

11 

H 

& 

Pi = P 

= m 

d( )/M = cL( )/a d0 = 

( )*/a 


= X = |[(h + t)t/ a 2] (E f /E c ) J 

The equilibrium equations (25) to (28), (31), and (34) reduce to 
the following equations: 


N' + Q + ap = 0 
Q* — N + aq = 0 


( 120 ) 


l^The effect of transverse shear stress deformation on homogeneous 
circular rings has been considered by L. Beskin in reference 6. 
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M' — aQ + am = 0 


og m = s - [M/(h + t)a] 


The stress-strain relations (%) , ( 59 ), (6 0), (6l), and (64) 
reduce to the following equations: 




Q = (h + t ) G c jj 3 + i(w* - v)J 


(1 + X)M = ± D*(P» + s/E c ) 

cl 


The load terms p, q, m, and s are given, according to 
equations ( 19 ) to (22), hy 


p - ( l + + C 1 - igrh 

« ■ ( x + + ( x - 
m ■ ^ [C 1 * ^y - ( x _ tt e ) pi ] 

(- ^>u - (l - 


1 

8 = 2 




'( 121 ) 


to ( 67 ) 


( 122 ) 


(123) 


(124) 


(125) 
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Ring sector acted upon by end 'bending moments .— As a first problem 
on circular rings, which illustrates the effect of transverse normal 
stress deformation, there is taken this "basic case for which, as is 
known, there must be the same stress distribution at all sections 
9 = Constant of the ring. 

According to equations (120) and (121), 


N = Q = 0 


M = ^ 




(126) 


= - ^A h + *)a 

J 


Equations (122) to (124) become 


v* + w = 0 
P + (w* — v)/a = 0 


(127) 


(1 + A)H> = D*p'/a 


(128) 


The significant result of this consideration is contained in equation (128), 
which may be written in the alternate form 


M = 


^ a 1 + \ 


B* 


1 (h + t)t E f 

2 a 2 E c 


V 


+ w 


(129) 


Thus, in this case of pure bending the transverse flexibility of 
the core is responsible for a reduction of the bending stiffness 

factor B* = ^ t(h + t) 2 Ef which is obtained exactly when E c = 0 

and practically when E c is of the same order of magnitude as Ef . 
Equation (129) shows that the reduction of B* is significant 
whenever E c is so small that the ratio E c /E f is of the same order 

of magnitude as the ratio (h + t)t/a 2 . 
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As a numerical example take the following values: h = 0.9 inch, 

t = 0.05 inch, a = 20 inches, and Ef/E c = 1000, for which 


1 (h + t)t Ef 1 0.95 X 0.05 

2 a 2 E c 2 400 


X 1000 = 0.0595 


indicating a reduction in tending stiffness of about 6 percent. 

Changing a from 20 inches to 10 inches changes the effect from 
6 percent to 2k percent. Changing Ef/E c from 1000 to 2000 increases 
the effect from 6 percent to 12 percent. Altogether it may "be said 
that this effect is of noticeable magnitude for some geometrically 
reasonable structures when the modulus ratio Ef/E c is of the 
order 1000 or more. Assuming aluminum face layers with Ef = 10? psi, 
this means that E c {g lO 1 *' psi, which is -well within the range of some 
present-day core-layer materials. 

Comparing equation (129) with the earlier formulas for the effect 
of transverse shear stress deformation, for instance with equation (ll6) 

in which a represents the plate radius and observing that G c ~ ^ E c , 

it is seen that the correction terms are of the same form, the difference 
being an appreciably larger numerical factor in the expression representing 
the shear effect. 

Closed circular ring acted upon by uniform radial load .— Having 
rotational symmetry, d/d0 = 0 and v = (3 = 0. Also set p = m = 0. 

The remaining equations permit the determination of the stresses in 
the face and core layers in a way which depends on the extent to which 
the load is applied to the outer (upper) and inner (lower) face membranes. 
Equation (12) becomes 


N = aq 


(130) 


From equations (121), it follows that 


^m = s ~ [M/(h + t)a] 


(131) 
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The stress-strain relations (122) to (124) give 


w 1_ 

a " C* 



(h + t)q _ aqf., X (h + t)C*"| aq 

12E c a C*|_ 3 l2E ca 2 J " C* 


(132) 


and 


M = 


D* s 
1 + X aE c 


(133) 


A closed circular ring subjected to a uniform radial load 
distribution q is stressed not only by a uniform axial force N = aq, 
as would be expected, but in addition is stressed by a uniform bending 
moment M, the magnitude of which is given by equation (133) . The 
explanation of this result is that for a ring with relatively soft 
core the circumferential stress distribution depends on the extent to 
which the external radial load is applied to the inner and outer forces, 
respectively. Roughly speaking, for a sufficiently flexible core layer 
the load q u goes predominantly into the outer face layer, while the 
load qj goes predominantly into the inner face layer. 

According to equations (46), in the present case for the stresses 
in the two face layers. 


S « - I " + iTTT 


M 


N7 = j N - T-i-r M 
1 2 h + t 


(134) 


According to equations (130) and (133) and in view of the definitions 
of D* and X, this may be written 

w _ a/ . 2Xs \ 

Nu -2{ q + TTTj 

Ni = —fq - - 2 - X s A 
1 2 V 1 + XJ 


( 135 ) 
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Combining next equations ( 131 ) and' ( 133 )* for the transverse 
normal stress in the core layer, the following expression is obtained: 


- 


1 + X 


(I36) 


For a specific example assume that the radial load is applied 
entirely to the inner face of the ring so that q u = 0 and, according to 

equation (125), 


» - C 1 - Wki 

6 = - s( x - 


' ( 137 ) 


With q and s given by equations (137), equations (135) and (136) 
tiecome 


-\k 

L h ± i 
2a 

;\ 

ji + x 



<^1 — 1 ^ 
H|CVJ 
II 

1 h + i 

2a 

(1 + 2 X)aq 2 
) 1 + X 

> 

( 138 ) 


“U - - ll* - J 2t- t )r7T 


It is seen that the flexibility of the core layer increases the clrcumfeir- 
ential stress in the loaded face layer in the ratio (l + 2 X) /(l + X) and 
decreases it in the unloaded face layer in the ratio l/(l + X), where X 
1 b defined by equation (119), compared with the equal values of these 
stresses when E c = 00. , 

Considering once more the numerical data under the section entitled 
"King sector acted upon by end bending moments," it is found, for 
instance, that the stress in the inner face layer may be about 6 or 12, 
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or 24 percent higher than the corresponding stress calculated without 
taking into account the transverse flexibility of the core layer. 

Ring sector acted upon by radial loads q u and q 1, uniform in 
circumferential direction and with vanishing resultant q .— Again it is 
assumed that d( ) /d0 = 0, m = p = 0 and now in addition that q = 0, 
so that, according to equation (125) , the only nonvanishing load 
term is s. Further, it is assumed that the ends 0 = ± a of the ring 
sector are free of stress, that is, W(ia) = Q(±a) = M(ln.) = 0. The 
ordinary theory of circular rings would then indicate the absence of 
deformations in the entire ring. In the present case there is found 
a type of deformation peculiar to the sandwich ring, which may perhaps 
be compared to the action of a Bourdon gage. 

Solving first equations (120) and (121) and satisfying the end 
conditions of the ring sector. 


N = Q = M = 0 


= 


= 8 




(l4o) 


The stress— strain relations (122) to (124) are then 


v 1 + w = 0 ] 


a(3 + w* — v = 0 > 


(l4l) 


3' = — s/Ec 


Assuming s independent of 0, from equation (l4l) there is obtained by 
integration, with constants of integration Aj_, Ag, and A^, 


3 = - Jj- 0 + Ai 

v = -a 0 + Aja. + Ag cos 0 + A3 sin 0 

= a -S- + Ao sin 0 — Ao cos 0 
E c ^ J 


( 142 ) 


V 
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As a specific example consider a complete ring, slitted radially 
at the section 6 = n, so that a = jt. Prescribe furthermore the 
symmetry conditions 0(0) = v(0) = w(0) = 0. Under these conditions 
there is obtained from equation (l42 ) 


F C P = — sS 


E c v = — as(0 — sin 0 ) 




( 1 ^ 3 ) 


E c w = as(l — cos 0) 


From equations (143), it follows that the radial slit, which is 
of zero width before the loads q u and qj are applied, opens und er 
the action of the loads to a width given by 


v(-n) - v(«) = 2na rg- = 2«a(l + h g + (l44) 


For a numerical example take a = 10 inches, h = 1 inch, t = 0.05 inch, 
E c = 10,000 psi, and q u = 20 psi, and obtain 


v(— n) — v(jr) = O.I32 inch (145) 


The foregoing three examples of ring analysis have been discussed 
in some detail, because they illustrate relatively simply the effect of 
transverse normal stress deformation in the theory of curved sandwich 
structures, without involving at the same time the effect of transverse 
shear stress deformation. 

Bending of semicircular ring by end shear forces .— Now a problem is 
considered in which both the values of E c and G c affect the result 
of the analysis. In the - equilibrium equations (120) and (121) all 
external load terms are set equal to zero and then, by integration and 
from the boundary conditions, that is, from 


N(± = M (± |) = 0 

Q(* I) = ^Qo 


( 146 ) 
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the following expressions for N, M, and 


Q are obtained: 


Q = Qo sin 0 ] 


N = Qq cos 0 


> 


M = — sQq cos 0 


(14-7) 


The stress— strain relations (122) to (124) become 


(1 + X./3)Qo cos 
Qq sin 9 = (h + 
— (l + X)Qoa cos 


9 = (C*/a)(v* + w) 
t)G c jp + (w* - v)/aj 


> 


9 = (D*/a)0' 


Integration of the last of equations (l48) gives 


D*0 = — a^(l + X)Qq sin 0 


(148) 


(1^9) 


where a constant of integration has been eliminated by means of the 
symmetry condition 0(0) = 0. Substituting equation (1^9) in the second 
of equations (l48). 



v) = 


Qo sin 0 
= Q Q sin 0 


t ) G c 


[jir+ 

D* |j + 2 


a^(l 


+ X) 


D* 

(h + t)t/Ef 


E r 


(150) 


Simultaneous solution of equation (150) and the first of 
equations (l48) for v and w gives as general expressions for v 
and v, 


v = A 0 cos 0 + A-l sin 0 + Aq cos 0 


w = A 0 sin 0 — (A^ + B)cos 0 + A 2 sin 0 


( 151 ) 
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where A]_ and A 2 are arbitrary constants of integration and A 
and B are found to he 


A = 


B = 


Qo a ' 


2D* 

Qoa3 


1 + 1 (h + t)t /E f + E f \ + p* 


V E C G Cy 


2D* 


^ p (h + t)t/Ef + Ef 


V E C G c 


a 2 C* 


D* 

a 2 C* 


( i+ i 

( i + 5 


(152) 


As further conditions, it is prescribed that v( 0 ) = vQQ = 0, which 
makes = A-j_ = 0 in equation ( 151 ). There remains 


v = A 6 cos 6 




w = A 6 sin G — B cos 6 


(153) 


Of particular interest are the values of w(n/2) and w(0), the first 
of these giving the radial deflection of the point of load application, 
the second giving the change of radius at right angles to the applied 
load. It is found that 



! (h + t)t/E f E f \ 

2 a 2 Ie c + G C J 




w( 0 ) 


= -B = - | 


Qo a 

D* 


1 

• 1 + 


]_ (h + t ) t /Ef 
2 a 2 \E C 



D* 

a 2 C* 



(15^) 


(155) 


Equations (154) and ( 155) contain the interesting result that, for 
this problem, transverse shear and transverse normal stress affect the 
outcome formally in nearly the same way. If the generally unimportant 
terms with D*/a 2 C* are omitted, which amounts to the usual assumption 
of circumferential inextensibility of the ring, then the effects of 
finite E c and G c occur in exactly the same way. 
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For a numerical example take h = 0.9 inch, t = 0.05 inch, 
a = 20 inches, E.p/E c = 1000, and Ep/G c = 2000. This giveB 


1 (h + t)t = 1 

2 a 2 16,800 

D* 1 (h + t) 2 _ 1 

a 2 C * l a 2 1770 

X = l (h + t)t % _ 1 

2 a 2 E c l6.8 

1 (h + t)t % _ 2 

2 a 2 G c 16.8 


The factors in brackets in equations (15^-) and (155) become 


1 + y , + y v + — — — (l + ^ ^ = 1.18 

16.8 16.8 1770 \ 3x16.8/ 

and 

1 + —^77 + —2— 1— (l + i— ) = 1.18 

16.8 16.8 1770V 3 x 16.8/ 


Thus, in the present example the flexibility of the core is responsible 
for an l8-percent increase of deflection-load ratio, and of this 12 percent 
is due to transverse shearing and 6 percent to transverse normal stress. 
Compared with these two effects the effect of circumferential extensibility 
of the composite ring is seen to be negligible. As a further numerical 
illustration, it is noted that reducing the ring radius a from 20 inches 
to 10 inches, with all other data unchanged, changes the l8-percent 
correction to a 72—percent correction. 

Bending of complete circular ring under action of two concentrated 
radial~f orces at 0 = ±k/ 2 .— The solution of this problem may be obtained 
by superposition of the solutions for the semicircular ring under the 
action of end shear forces Qq (equations (l46) to (155)) and under the 
action of end loading moments Mq (equations (126) to (129)). 
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The first step consists in determining. Mq in terms of Q q such 
that the sum of the 0’s from equations (129) and (1^9) assumes the 
value zero for 9 = it/2, that is, the value of the superimposed tending 
moment at 9 = it/2 must make the tangent to the deflected ring at 
this point horizontal. Combining equations (129) and (1^9) in this 
manner , there is obtained 


it 1 + 1 c, 1 + 1 „2 Q _ /-v 

2 d* aM 0 p* a Wq - U 


or 
\ ' 

M 0 = (2/*)aQo (156) 


It may be noted that equation (156) is a further case of a statically 
indeterminate problem where transverse shear and normal stress flexibility 
do not affect the internal force and moment distribution but affect only 
the state of deformation of the structure. 

Further , the radial deflections w(«/2) and w(0) due to the 
action of Mq are calculated, in order to combine them with equations (152) 
and (153). Integrating equations (129) and (127) with the boundary 
conditions v(0) = v(jt/2) = 0, there is obtained for the displacements 
due to 


D*w = —(l + l)MQa^l — ^ cos 0^ 
D*v = (1 + l)MQa 2 (e - £ sin 0) 


(157) 


and, in particular. 


D*w('0) = (1 + l)Moa 2 (J - l) 
D*wg) = -(1 + l)Moa 2 


( 158 ) 
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Combining equations (158) with equations (154) and (155) and taking Mq 
from equation (156), there follows for the resultant displacements 



. , (h + t)t Ef 

where kn = 5 tv has been put as a further abbreviation. 

2a 

Equations (159) may be written in the alternate form 


v(*/2) = 0.149 \ 


1 + \ + 5.29Uc + 




(160) 


w(0) = -0.137 


Qoa3 


1 .+ X + 3 . 65 fx 

Q. - -^-(l + W 1 

L 

U a 2 C* ' 3/J 


(161) 


When X = Xg = 0 and when the composite ring is assumed axially 
inextensible, which amounts to putting D*/a 2 C* = 0 inequations (l60) 
and (161), then equations (l60) and (l6l) reduce to well-known results 
of circular-ring analysis. 

Comparing equations (160) and (l6l) for the closed circular ring 
with equations (154) and (155) for the open semicircular ring, it is 
noteworthy that for the semicircular ring X and Xq occur with equal 
weight, while for the closed circular ring the influence of Xq is 
considerably greater than the influence of X. Thus, for the closed 
circular ring the effect of transverse shear deformation is much more 
important than the effect of transverse normal stress deformation, while 
for the open semicircular ring both effects occur in a much more nearly 
equally important way. 

For a numerical example of the use of equations (l6o) and (l6l) 
take again the values for the numerical example given in the section 
entitled "Bending of semicircular ring by end Bhear forces." This 
gives for the expressions in braces 
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and 


1 . + g x 5.29 ^2£ 

16.8 16.8 1770 


1.69 


1 + 


1 . 2 X 3.6^ 

16.8 16.8 


+ 1^5 
1770 


1.50 


Thus, while the effect of transverse stress deformation for the open 
circular ring amounted to 18 percent, the corresponding corrections 
for the closed ring are 69 and 50 percent, respectively. 

The next step in the analysis of sandwich— type circular rings 
would be the general, solution of the system of equations (120) 
to (124) for arbitrary load distributions. This, evidently, is 
possible and further specific examples of interest might be analyzed 
on the basis of the general solution. Such extension of the work of 
this section is, however, left for future considerations. 


Circular Cylindrical Shells 

In this section the general system of equations of part I of this 
report is restricted to the equations of the theory of circular cylindrical 
shells. The treatment of sandwich— type shells of this kind is shown to 
be not appreciably more difficult than the analysis without the effect 
of transverse shear and normal stress. 

As specific examples some problems of rotationally symmetric 
deformations are treated. In particular the influence coefficients 
are obtained for a semi— infinite shell acted upon by bending moments 
and transverse forces at one end of the semi— inf inite shell. With these 
influence coefficients an explicit solution is obtained for the problem 
of the infinite circular cylindrical shell acted upon by a pressure 
band of zero width . 
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In the general equations of the problem there are set for the 
relevant coordinates and variables , 


*1 = 

ae 

= x 

= 

r 2 = 

00 

H 

II 

c? 

N 22 = N x 

Ql = 

00 

II 

II 

1— 1 

M 22 = 

- Mr 

U]_ = V 

d 

11 

CVJ 

d 

Pi = 

P0 

P 2 = Px 

n 

= 

“x 

►d 

H 

II 

►d 

CD 

a 

11 

OJ 



£ 

II 

■J 8 

II 

H 




N 12 = n 21 = N x0, 




M 12 = ^1 = Mx0 



(162) 


The equilibrium differential equations (25) to (28), (31), and (34) 
become 


dN x X 8N x e 

Bx a d0 


+ p x = 0 


^ N xe 1 

dx a SO a 


+ 


Pe 


0 > 


52s 

dx 


+ i 

a S9 


SQe N ( 


- -r + q = 



(163) 


1 ^c0 
dx a S0 


dMge ,1^9 
dx a bO 


Q x + — 0 

% + Dig = 0 


( 164 ) 
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°u = s - 


Me 


(h + t)a 


(165) 


The stress-strain relations (5 6), (59), (60), (6l), and .(64) to (67) 
become, with X2 = ^12 = 0, = i[(h + t)t/a 2 J(Ef /E c ) = X, 


(1 ♦ i)* e - v* = c[i g ♦ I ♦ I^la] 

"z - VN e * °*(i) 

2(1 + v)h i9 . 0^ ♦ i 


(166) 


Qe = (h + t)G c ^ e + 1 
Qj = (h + tjGc^ + 


(1 b \ 

\a Se~ aE c y 


(1 + X)J% — VMg. = D* 


Mj — VMq = D*^ 


2(1 ♦ v)H* - d*(| 2 + 1 |i) 


(167) 


(168) 


When G c = E c = ~ (and therewith X = 0) equations (163), (l64), 
(l66), (167), and (l68) reduce to the known system of equations in which 
deformations due to transverse stresses are neglected. The solution of 
the present system of equations is not essentially more difficult 
the solution of the system with G c = E c = «o. in particular also here 
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there may he obtained a trigonometric double-series solution, as a 
generalization of Wavier' s solution for the flat plate (references 7 
and 8) . 

For this trigonometric double-series solution there is set. 


q = ^ qnm ain me sin nx/Z 
P 0 = ^ ^ P0nm 008 me sin n*/ 1 

Pi = ^ Prrrm sin m B cos m/Z 

%: = ^ gyrrm a in m 9 cos nx/Z 

“0 = ^ ^ ®0mn C0B me ® in nx/z 

s = ^ Smn sin me sin nx/Z 


(169) 


w = VTrm sin me sin nx/Z 

v = Virm cos me sin nx/z 

-II Ujga sin me cos nx/Z 

’-II Pxmn fl ln COB nx/Z 

Pe = ^ ^ Pemn COB m0 Bl -n nx/z 


( 170 ) 
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Sc = 
Qe = 


II 

II 


Q^jjh sin m0 cos nx/Z 
Qflmn cos me sin nx/Z 



( N x * ~ ^ ^ ( W xmn^ %mn) a ^ n s ^ n n^/Z 

N x0mn 008 me cos 1X1 A 

M e) = ^ ^ (Mymn. %nn) 8in m0 sin “A 
M^g = ^ ^ ^©nm cos m© cos nx/Z 


(171) 


When equations (169) to (l7l) are substituted in equations (163) 
to (l68) there remains for every value of m and n a system of 
13 simultaneous equations for the 13 Fourier coefficients which occur 
in equation® (170) and (171). 


A system of only five simultaneous equations for the five Fourier 
coefficients in equation (170) is obtained if first equations (163) 
and (l64) are reduced to five equations for the five unknowns w, v, u, 
3 X , and by means of equations (l 66) to (l68). 


For the present, the task is not carried out of obtaining the 
deformation and internal stress Fourier coefficients of equations (170) 
and (171) in terms of the Fourier coefficients of the load terms in 
equation (169). Instead, the axisymmetrical case, to which 
equations (169) to (l7l) reduce when sin m0 and cos m0 are 
interchanged throughout, and then only the terms for m = 0 are taken, 
is treated separately. 

Axisymmetrical deformation of circular cylindrical shell .- In 
equations (163) to (l68) set 


d( )/Se = 0 

a( )/dx = ( )* 

W X0 = % = %X = 0 >■ 

V = 3 e =■ 0 

m e = p© = 0 


(172) 
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and then the following system of equations has to he dealt with: 


N i’ + Px = 0 

Ox' - (%/a) + q = 0 


(173) 


Mx* - Q± + % = 0 

°£m = s - %/(h + t)a 


(174) 



(175) 


Qj = (h + t)G c (P x + w’) 


(176) 


(1 + X)% - VMjj. = D*s/aE c 
Mj - V% « D*p 2 ' 


(177) 


The system of equations (173) to (177) may he reduced to two 
simultaneous equations for 0 2 and Q^., as follows: First, 

express Mjj. in terms of 0 X hy means of equation (175) and substitute 

the result in equation (174). From the first of equations (177) it 
follows that * 


% = 


V 

1 .+ X 


Mjc '+ 


D*s 


(1 + X)aE c 


(178) 
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and this, introduced into the second of equations (177), gives 




X)D* 


1 + X - V 2 


VD»s 


aE c (l + X - v 2 ) 


(179) 


Equation (179) is introduced into the first of equations (174) and, 
restricting attention to shells of uniform section properties, there 
is obtained 


/ 


(1 ± _ X ) D1_ p 
1 + X - V 2 


~ Ox = - 


VD*S 

aE c (l + X - v 2 ) 


(180)* 


To obtain the second of these equations, first, introduce into 
equation (176) the value of w’ which follows from equation (175), 
giving 


gx 

(h + t)G c 



" VN X * - 


C*(h + t)q* ” 
- 12aE c 


(.181) 


In equation (l8l), N 0 * and N x ' are taken from equation (173) and, 
after slight transformations, there is obtained 



0 * 

(h + t)G c 


+ 3 X 



(182)* 


Comparing equations (180) and (182) with the corresponding 
equations without the effect of transverse shear and normal stress 
deformation, it is seen that the effect of transverse normal stress, 
which is represented by X, merely somewhat modifies some of the 
coefficients of the left sides of the corresponding system of equations 
with E c = co. In contrast to this, the effect of finite G c is to 

introduce a new term into the left sides of these equations. ThiB new 
term may be of appreciable importance, as will be shown. 

Having solved equations (l80) and (182), and Mg are obtained 
from equations (179) and (178), respectively; Np follows from 
equation (173) In the form 


N 0 = a(Qj t ’ + q) 


(183) 
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and v follows from equation (175) In the form 


w = (a/C*) (1 + 1/3^0^ + aq + v 




(184) 


The following examples illustrate the use of equations (178) to (l84) . 

Infinite circular cylindrical shell with periodic load Aistrlhution .- 
In specialization of equations (169) to (171), set 


q = q^ sin MX 
8=8^ sin px 
ftc = Pxm COB MX 
= %m 008 ^ 


w = W M si - n MX 
u = cos mx 

Px - Pxm cos mx 




= Sim cob 

®x = ®xm s ^ n MX 

N 0 = N 0 4 sin MX * 
He = m xm sin MX 

' = ^m 8in M* 


(185) 


(186) 


t 


(187) 
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By introducing equations (185) to (187) into equations (l80) 
and (182), two simultaneous equations are obtained for the ampli— 
tudes and 0 X(i , as follows: 


(l + X)D* 2/3 n VD* 

1 + X - v 2 4 ^ + ^ + 1 + X - v2 aE 


c 


c* 


H) 


M 2 + 


2tE, 


a 2 (h + t)G r 




(188) 


Q *H = - + vPxn/a) 


To simplify the further discussion, by setting in equation (l88) 
m^ = s^ = p X(i = 0, there is obtained for P X)a and 


1 + X - v 2 % K 


(1 + X)D* u 3 




(189) 


The quantity K is given by 


K = 


n X 2 
1 + 3 + ^ 


Z 2 t 


E. 


f 1 + X - v‘ 
+ 




~i-l 


it 2 (h + t)a 2 G c (1 + x) iA (h + t) 2 a^| 


(190) 


where use has been made of the relation i_i = n/l. In equation (190) 
the term X/3 will usually be of little importance. The other two 
variable terms represent the effect of transverse shear deformation 
and of shell curvature, respectively. When the radius a is so large 

that Z^/(h + t) 2 a 2 « 1, the shell behaves under the action of the 
given load essentially as a plate strip. The effect of transverse 
shear is important as soon as the term (2/n 2 )(2 2 /a2)(t/(h + t))(Ef/G c ) 
is not small compared with 1. 

Before evaluating a numerical example the following further 
formulas which are readily obtained from equations (179 ) , ( 183 ) , 
and (l8U) are listed: 
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M>cn = 


N 0H = a^ 1 ~ K ) 

t 1 - ( x + D 


? 

a 4.4 
w h = — 


( 191 ) 


Equations (l9l) show that in this problem not only is the deflection 
increased because of the effect of transverse shear, and with that 
the hoop stress resultant Ng^, but now also an effect is found on 

the bending-moment distribution M^, in the opposite sense. The effect 
of transverse shear is to reduce the magnitude of the bending moments 
in the shell. This result is in contrast to what was found for the 
examples which were worked out in the sections on plate analysis and 
circular ring analysis and iB therefore of particular significance. 

Equation (191) for w^ may be compared with the corresponding 
expression for a simply supported plate strip of width Z, with 
sinusoidal load. The result for this case must follow from equation (191) 
in the limit a^-oo and agree with equation (96), which was previously 
obtained. To compare the last of equations (191) with equation (96), 
the last of equations (191) is written in the form 


w,, = 


1 + X - y 2 n\ k q + /JL ) 2 1 ± x 1 ( h + t)t Ef 

(1 + X)P* W W 1 + x - v 2 2 Z 2 . g c 




(1 + X)D* 


1 + 


X 

3 


0 


2t 




h + t G 


1 + X — v 2 


% (192) 


Equation (192) reduces to the equivalent of equation (9 6) if in 
it a — ^ 00* 

From a comparison of equations (192) and (96), it is further 
concluded that the correction due to transverse shear is greatest 
in this case when a = <», so that, in this case, the curvature of' 
the shell tends to reduce the additional shear deformation, below 
the value obtained for the simply supported plate strip. 
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For a numerical example firBt take h = 1 inch, t = 0.05 inch, 
a = 10 inches, Z = 20 inches, Ef/G c = 200, Ef/E c = 100, V = 1/3^ 

and X = 1 100 = °-°25. The factor K of equation (190) 

he comes 


K 


1 


1 + 0.008 


2 400 X 0.05 „ . 4(1 - 0.09) 160,000 1 - 1 

100 X 1.05 J ioo X 1.1 


= (1 + 0.008 + 3.86 + 5^.5) 


= 0.01685 


while without transverse shear and normal stress deformation 


( k ) g c = E c = o . = (1 + 5 ^. 5) -1 = 0.0180 

The correction in tills case amounts to about 6 percent. 

Changing the moduli ratio to Ef/G c = 2000, Ef/E c = 1000, 


K = (1 + 0.08 + 38.6 + 54.5) = 0.0106 


instead of K = O.OI685. The correction in this case amounts to 

— T ' ^Q^QiQg > | Q1Q ^ x 100 * 70 percent . Thus again there is a case where 

omission of the effect of transverse shear deformation would give 
results which could not he used. However, it is noted that the effect 
of transverse normal stress deformation is quite nma! 1 »,nd may here 
safely he neglected. 

If the foregoing values of K are introduced into equations (191), 
it is seen that the percentage corrections apply to the bending-moment 
value directly hut that for hoop tension and radial deflection the 
corrections are very small indeed. In fact, in order that there he 
appreciable corrections due to transverse shear on hoop tension and 
radial deflection, it is necessary that the half wave length of the 
sinusoidal load q he so small that K is at least of magnitude 0.25 
or more. 


/ 
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A case of approximately this kind is obtained if the half wave 
length 2 is changed from 20 inches to 10 inches and the moduli ratios 
are again taken as Ef/E c = 100, Ef/G c = 200. Then, 


K = (1 +0.008 + O.965 + 3.4l) -1 = 0.1865 


whereas 


(K) Ec=Gc=<a = (1 + 3 ad- 1 = 0.227 

The percentage change of K and therewith of Mj is slightly more 
than 19. The percentage change of Ng and w is about 

The foregoing numerical examples show that the effect of transverse 
shear may be significant in cylindrical sandwich-shell analysis and that 
moreover its magnitude will not in general be predictable by the a nal ysis 
of an equivalent flat— plate or straight— beam problem. 

For the infinite circular cylindrical shell with 
load q = q 4 cos px the essential results are given by equations (190) 

and (191). These results may be extended directly to the loading 
condition 


q = ^ q n cos p n x 

t 

Mix = n *A 




(193) 


superposition, from equation (191) the following formulas are 
obtained: 


% - Y (in/Mn 2 )^ cos p n x 


% = a 2^ TnO- ~ Kq) 008 M n x 

v = (a2/C*) Y_ qn[l ~ (1 + V3)K n ] 




(19M 


cos 
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The values of are obtained from the formula 


K n = 


1 + ^ + 


z 2 t % . 4(i + x - v2) 

. n ' 11 


3 n 2 * 2 (h + t)a 2 G c (1 + x) Ai^ (h 


1-1 


t) 2 a s 


( 195 ) 


Having the solution for the infinite shell with periodic load 
distribution, it will be only necessary to add to this the general 
solution of the differential equations without external load terms, 
in order to obtain the complete solution for any edge condition of 
the axisymmetrically stressed circular cylindrical shell of finite length. 
This additional solution will now be obtained. 

Finite circular cylindrical shell acted upon by edge moments and 
forces .— To solve equations (l80) and (182) with right-hand sides equal 
to zero, equation (182) is differentiated twice and p x * ’ ia 
substituted from equation (l8o). This gives 



1 

(h + t)G c 


Q, 


• t 


+ 1 + 1 - V 2 
(1 + X)D* 


Qx - 0 


or 


Qy 17 - 2m 1 2 Q x " + Umg 4 ^ = 0 


(196) 


where 



C* 


+ X/3)(h + t)G c 




(1 + X/ 3 )(h + t) 


!f 

&C 


(197a) 


**-N| 


I c* 


1 + X - V e 


a 2 D* (1 + X)(l + X/ 3 ) 



1 + X - V 2 
+ X)(l + X/ 3 ) 


( 197 b) 
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The auxiliary equation corresponding to equation (19 6) is 


r.^ — 2m^r^ + = 0 


(198a) 


or 


r 2 = m3 2 ± m-^ - 


(I98h) 


The solution of equation (196) occurs in two different forms, depending 
on whether r 2 of equation (198b) is real or not. According to 
equations (197) and (198b), r 2 is complex as long as 


m^ < kmq 1 


or 


(199a) 


Jl_[ tEf * 

a 4 _(h + t)(l + X/3)G C 


2 


< 


4 [l - v 2 /(l + 1)1 

(h + t) 2 a 2 (l + X./3) 


To clarify this condition, neglect X (which is of very little importance 
here) and equation (199a) then becomes 


?£< 


2 & 

t 


When equation (199) holds, a quantity k may be defined by 


(199b) 


k = 




( 200 ) 


and the four roots of the characteristic equation are k, k, — k, 

and -k, where a bar indicates the taking of conjugates. The solution 
of equation (196) may be written ’ ' 

Qx = Cqe ^ + C^e ^ + 026^ + C2e' 


( 201 ) 



68 


NACA TR No. 1832 


Where equation (199) does not hold, which is the case for very small 
values of G c /Ef only, all four roots of equation (198a) are real and of 

the form 


k 1 = \j + \j 

& 

II 

k 

H 


k 3 = Sj*! 2 - ' 

yjm^ — 4m g 2 

k 4 = -k3 



( 202 ) 


gnfl the solution of equation (196) can he taken in the form 


kix — k-ix koX — kox 

= A-je + Age + A^e J + Aj^e J 


(203) 


Before applying either solution to a specific problem, there are 
noted the following relations which follow from equation (200) 


kk = |k| 2 = 2mg 2 
k + k = \[2^ \j ta.^ 2 + 2ni2 2 


(204) 


Semi— infinite shell acted upon by edge bending moment and shear 
force .— There are the following boundary conditions. 


Mjc(0) = (1 - + ^)B g- P x * ( 0) = M o 
1 + X - V 2 

Q x (0) = Qo 

while for x = <» these same quantities vanish. ^5 


(205) 


^5For the same problem without the effect Of transverse shear and 
normal stress, see reference 9. 
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Of particular interest in this solution are the values of 
deflection w(0) and change of slope 0 X (O) at the section where 

the loads Mq and Qq are applied. 16 

Taking first the case Ef/G c < 2a/t for which equation (201) 
applies, it is seen that the conditions at infinity require that 


C 2 = C 2 = 0 


( 206 ) 


so that 


— kx _ — kx 


(207) 


The values of may he obtained "by integration from equation (l80) 

in the form 


(1 + X)D* 

1 + X - V 2 


3x 


e -kx ^1 — kx 

* V 


(208) 


where two constants of integration have “been discarded to satisfy again 
the conditions at infinity. 

With equations (207) and (208) there is obtained from the "boundary 
conditions (equations (205)) that 


Ci + Ci = Qo 


(Ci/k) + (Cj/k) = -Mq 


(209) 


^Without transverse shear and normal stress deformation these 
relations are 


w(0> - 


Mq +' 


4a 2 D-» 


(l - v 2 )c* 


Qo 


Mo) 


= \[^ 
\l D*i 


D*C* 


A M (l - V 2 )c* „ 

Qo + 2 \ — — o *^0 

M 4a 2 D* 


in agreement with equations (236) of reference 9, where the homogeneous 
shell is considered. 
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This determines and in the form 


c kQp + |k| Mp 
1 k - k 


- = kOp + |k| 2 Mp 

1 . k - k 


( 210 ) 


Equation (210) is introduced into equation (208) and there is 
obtained as the first of two "influence coefficient" formulas 


(Vx X - P v g P:t(0> " - + O' + k)Mo] (211) 

The second of these formulas follows from equations (l84), (207) 
and (210) in the form 

1 r V 3 ^ = -1 h 2 m o + (k + *4 (2i2) 

Equations (21l) and (212) ma^ he written in more explicit form, 
using equations (204) and (197) . The results are 


M0) = - 


\/c*D* 




(1 + x/3)(l + x - V 2 ) 


1 + X 


Qo 


+ 



1 + X - V £ 


a 2 D* (1 + X)(l + X/: 


+ (C*/a 2 G c ) 

(h + t)(l + X/3) 



( 213 ) 



NACA TN No. 1832 


71 


and 



Neglecting the generally small effect of finite E c in 

equations (213) and (214), that is, putting 1=0 in these equations, 
there may he written instead 



Equations (215) and (216) contain the noteworthy fact that the 
correction factors for the effect of transverse shear are independent 
of the ratio t/h of face— layer thickness to core thickness. The 
complete formulas of course must and do contain the influence of the 
core thickness h. 

It is further noted that, while equations (21l) to (2l6) have been 
derived for the case that m^ < 4nirA, for which the complex solution 

holds, they are also valid, as is readily shown, when > m^. 

Comparing -equations (213) and (214), and (215) and (216) with the 
equations listed in footnote 1 6 it is seen that: (l) The effect of 

transverse shear modifies the deflection due to Qq and the rotation 

due to Mq hut not the other two coefficients, (2) the effect of 
transverse normal stress enters all four coefficients hut only in a 
minor way, and (3) the reciprocity relation that the deflection due 
to Mo is the same as the rotation due to Qq is carried over from 
the theory without the extra effects. 
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For a numerical example the following data are chosen: t = 0.1 inch, 

h = 1 inch, a = 10 inches, Ef/E c = 100, Ef/G c = 200, V = 1/3, 

and X = ^ 100 = 0.055, and, from equation (197), 


mi 


_1_ 0.1 X 200~ 

10 \j 1.018 X 1.1 


0.426 


“2 " 


1- QfQ9 

100 X 1.21 X 1.018 


0.294 


Then, according to equation (204), 


|k| 2 = 0.173 


k + k= 0.182 + 0.173 = 0.84 


while without transverse shear deformation (mi = 0) the value 

of k + k = O.59. According to equations (211) and (212), the effect 
of transverse shear in this case is to increase the rotation due to 
the edge moment in the ratio 0.84/0.59 = 1.42, an effect of 42 percent. 
The same increase is found for the deflection due to the edge shear 
force. Rotation due to the shear force and deflection due to the 
moments are practically unchanged. Likewise, the effect of transverse 
normal stress in this case is of negligible importance. 

As a further numerical example there is chosen t = 0.05 inch, 
h = 1 inch, a = 20 inches, Ef/E c = 1000, Ef/G c = 2000, 

and X = i 1-^O.Cg iqqo = O.O 65, and, from equation (197), 


m l 



0.05 

1.022 X 1.05 


2000 


0.483 


1 n 4| 1-0.0 9 

\J 1.05 X 20 ^ 1,022 


0.218 


m2 = 
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From equation (204) then 


|k| 2 = O.O95 

k + k= NfFxf 0.234 + 0.095 = 0.82 


while without transverse shear deformation (mj_ = 0) the value 

of k + k = 0.44. Thus the effect in this case is to increase edge 
rotation due to edge moment and edge deflection due to edge shear 
force in the ratio 0.82/0.44 = I.87, an effect of 87 percent. 

Infinite circular cylindrical shell acted upon hy transverse line 
load .— Calculation is restricted to the determination of deflection and 
bending moment at the section x = 0 where the line load of intensity 2Qo 
is assumed to act. The result of the foregoing paragraph may be used as 
follows. Consider the infinite shell cut in two parts at the section x = 0 
and assume a bending moment Mq of such magnitude that the slope 0 X (O) 

is zero. According to equation (211), this gives 


Mq = - 


Qq 


Qo 


k + k \J 2m x ^ + 4m2 2 


(217) 


and > therewith 


, (k + k) 2 - |k| 2 _ g mi 2 + gm, 2 

(1 + x/3)a 2 . k + k \jaii! 2 + kmz 2 


(218) 


Equations (217) and (218) become, with equations (204) anti (197), 



% 

\l C * 1 + 1 - V 2 , C^/a 2 1 

\ N/a^D* (l + X)(l + X/3) (h + t)G(> 1 + X./3 


( 219 ) 
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w(0) 


(1 ♦ l> 2 \ 

C* 1 + X - v 2 ! C*/a 2 1 

/ a 2 D* (1 + X)(l + X/3) ' (4 + t)G c 1 + X/3 

’ i 

2 \ 

C* 1 + X - v2 ( c*/a 2 1 

a 2 D* (1 + X)(l + X/3) ' (h + t)G c 1 + 4/3 


( 220 ) 


To give these formulas a somewhat less unwieldy appearance, the 
effect of finite E c , that is, X = 0, may again he neglected i as is 
permissible in most cases j and there may he written 



( 221 ) * 


n It ®f 

, n , - - N&Iv! a? N/l - v2 5 ^ % 

k t \/(h + t)a | l t ff f 
\| 2 \/ 1 - v 2 & G ° 

Some numerical examples are given as follows. 

Taking t = 0.1 inch, a = 10 inches, Ef/G c = 200, and V = l/3, 
transverse shear deformation reduces Mq to l/\/2.05 times the 
value which holds when G c = <», that is, there is about a 30 percent 
reduction in Mq. At the same time the deflection under the line 

load is 3.05/ \j2.O5 = 2.14 times what it is when G c = »; that is, 
there is an increase of about 115 percent in w(0). 

Taking t = 0.05 inch, a = 20 inches, and Ef/G c = 200, , Mq is 
decreased by a factor \j 4/5 = O.89, while w(0) is increased by a 
factor 1.5/ \jl.25 = 1.34. 

Taking t = 0.05 inch, a = 20 inches, and E f /G c = 2000, Mq is 

decreased by a factor l/ \] 3.62 = 0.526, while w(0) is increased by 
a factor 6.25 / \/ 3 - ^2 = 3*29. 


( 222 )* 
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Equation (220) for w(0) may "be compared with equation (ll6) for 
the circular plate of radius a. This comparison shows that, while for 
the plate "both the ratios t/a and (h + t)/a enter into the correction 
factor, the correction factor for the cylindrical shell contains the 
ratio t/a only; that is, the corrections (hut not the results) are 
independent of the ratio of face— layer thickness to core thickness in 
this case of a cylindrical shell. 


Spherical Shells 

In conformity with customary usage, the following notation is 
introduced : 


= a 0 

| 2 = 

= 1 

ctg = sin 0 

bd 

H 

tl 

II 

P 

%1 = N0 

E 22 = N 0 

II 

H 

II 

$ 

% = Q0 

3 > 

II 

& 

H 

II 

^22 = % 

M i2 = ttye 

H 

II 

I >2 = P& 

h 0 

II 

ID^ ~ 

u^ = u 

u 2 = v 

Pi = P0 

P 2 = Pe 


(223) 


Attention is here restricted to problems with rotational symmetry 
and the following relations are used: 


d( ) fie = 0 

N 0e = % = = Pe = ^ = v = 3 0 = 


•1 


( 224 ) 


The differential equations of equilibrium (25) to (28), (31) 
and (34) become, setting 3 "" ' 3 


d( )/&0 = d( )/d0 = ( )• 
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(sin 0N0^' — cos 0N 0 + sin 0Q^ + a sin0 = 0 "(225) 

(sin 0Q^’ - sin 0^N^ + N 0 ^ + a sin0 q = 0 (226) 


(sin 0M^y — cos 01^ — a sin 0Q^ + a sin 0 = 0 (227) 


+ (^0 + + t)a - s = 0 


(228) 


The stress— strain relations (.56), (59) to (6l), and (64) to (67) 
■become, if there is set in accordance vith equation (63) 


. 1 (h + t)t E f 

~ K 22 ~ 2 a 2 E c “ * 

(229) 

- !> ■ °*( u ' : w * Ma* i) 

(230) 

- i)v - <K na °v Uv + £) 

(231) 

(H + tjG,.^ +»' -«) 

(232) 

(- (v-x)Me -¥(^’ '+£) 

(233) 

(V - X)M0 > cot 0 + £) 

(234) 
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There is first given a simple special solution of this system of 
equations and then a generalization is obtained of the two simultaneous 
equations for Q0 and 00 which are fundamental in the theory of 

homogeneous isotropic shells. 

Uniform stress distribution in a spherical shell .— Set in 
equations (225) to (23*0 p0 = m0 = 0 and assume that Q^ } M^ } 

and Mg are independent of <f>. From equation (225) it follows that: 


N0 = N 0 = N 0 

Q0 = 0 


(235) 


From equation (226) it follows then that 

N 0 = | 


(236) 


and from equation (227) it follows that 


M0 = Mg = Mq 


(237) 


Equation (228) gives 


CT^m = 8 - 2M0/(h + t)a 


(238) 


In equations (230) and (231) set u = 0 for reasons of symmetry and 
obtain 


(1 + | X - v)n 0 » e w + 


t(h + t) Ef 


6a 


Er 


(239) 
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s 


or, with Nq from equation (236) and X from equation (229), 


(C*/a)v = |(l - V)aq (240) 

Equation (232) is identically satisfied when 00 = 0. Equations (233) 
and (234), in conjunction with equation (237 ) , give 


(1 + 2X - V)Mq = (D*/a) (s/E c ) = | s 

c 


or 


Mq = 


(h + t)aX 

1 + 2X - v 


(241) 


Then, from equation (238), 




Jl ~ v )s 

1 + 2X - v 


(242) 


Equation (242) may "be compared with equation (136) for the circular ring. 

According to equation (46), there are obtained from equations (236) 
and (24l) the following expressions for the stress resultants in the 
outer ("upper") and inner ("lower") face layers: 




+ - ^ 

1 + 2X 



y 



Is 

1 + 2X - V 



(243) 


Comparison of these results with the corresponding results for the 
circular ring (equations (135)) shows that for given values of q 
and s there is a greater difference "between N u and N^ in the 
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spherical shell than there is in the circular ring, the reason "being 
the relatively larger influence of the s-term in equation (243) . 

For a specific example, it is again assumed that the radial load is 
applied entirely to the inner face so that q u = 0 and, according to 

equations (20) and (22), 



Substitution of equation (244) in equation (243) gives 





h + t \|j _ /-> h + tV~ aq z 
2a A " v 1 ~ j ~r 

= fl - h + t\ ai n 1 + 4X - V 
\ 2a J 4 1 + 21 - V 


1 - V 

1 + 2X - V 




J 


(245) 


As a numerical example, taking X — 0.0595^ us in the example given in 
the section entitled "Closed circular ring acted upon by un i form radial load 
and V = 1/3, it is found that the factor in N z which contains the effect 
of the core flexibility is (l + 0.36)/(l + 0.18) = 1.15. Thus, where 
for the circular ring there was a. 6-percent stress increase, there now 
is a 15— percent stress increase. 

Reduction of axis^iutuetrical problem to two simultaneous equations 
for Q0 and Pri. — The fundamental results of reference 10 for 

homogeneous shells may be readily extended to sandwich shells as 
follows: 

Equations (225) and (22 6) are used to express N^ and Ng in 
terms of Q0. 

w 0 = cot 0Q0 + F x (0) (2 46) 


(247) 


% = 00* + P 2 (0) 
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In equations (246) and (247) the functions F^ and F 2 are given "by 

Fl = ^ J* ^q cos 0 - p0 sin 0)sin 0 d0 (248) 


F 2 = " CO g 0 [( aln 0 F l)' + ap 0 8in $\ (249) 

Next the displacement components u and w are expressed in terms 
of 00, by means of equations (230), (231), (246), and (247)-. 

Subtraction of equation (231) from equation (23O) gives 


£(u*. “ cot 0 u ) = (1 + V)(N0 - N 0 ) 


= (1 + V)[-(Q0' - cot 000) +. F-l - Fgj (250) 


Equation (112) is integrated to 


(C*/a)u = -(1 + v)(Q0 + F^) (251) 


where ia given hy 


^ J sin 0 


(252) 


Equations (25l) and (252) are introduced into equation (231) and the 
following expression is obtained for w: 


(C*/a)w = (1 + 1/3) (cot 0Q0 + Q0 f ) + 


(253) 
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where F 4 is given "by 


F k = - ^ ^ + (1 + v)cot 0 F 3 + (l + ^)F 2 - (v - ^F X (254) 


Equations (251) and (253) are introduced into equation (232) for 
and the first of the two simultaneous equations for and 30 is 

obtained in the form 


(h + t)G c 


C* 


(l + ^ (cot 0 Q0 + Q0*)‘ + F 4 » 


+ (1 + V) (Q0 + F3) 


which may he rearranged to read^-7 


cot 0 


— cot^0 — — — 
! + 


~ 1/3 , Cf 

/ 1 _ 


X/3 G c (h + t)(l 


(1 + X/3)] ( 


T7TR ^ = " 5(0) 


(255) 


The function F^ is given hy 


d + v ) f 3 + F 4 

5 ’ — " l/i ' ' 


(256) 


Introducing the operator 


L = ( )' + cot 0( )' - cot 2 0( ) 


ITWhen X = 0 and Gc = » and when no external loads are present, 
this equation checks with the first of equations (g) on page 469 of 
reference 9- 
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equation (255) may finally also "be written 


(L - U1 )Q0 + YTTfi H * r 5 W . (257)* 

where 


^1 


(_ 2t 5c 

h + t G c 


_ v + JL_^ 

3/1 + X/3 


(258) 


The second of the two simultaneous equations is obtained somewhat 
more directly as follows: Write equations (233) and (23*0 in the form 


** = r~^T^oT V) [ (1 + x) ^' + (v " x)o,t ^ H 

+ (1 + V)s/E c ] ' (259) 

% Rl + *)cot 0 30 +(v - 1 ) 30 ' 

1 - V 2 + 2X(l + V) L F ^ 

+ (1 + v)s/e c J (260) 


Introduce equations (259) and (260) into the moment equilibrium 
equation (228) and obtain 


D*/a 


1 - V 2 + 21(1 + v) 


+ cot 0 30* - (cot 2 0 + j — ^ 30 + (1 + v)s'/E c 

+ (1 + v)s*/E c J - aQ^ + am^. = 0 (261) 
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Again, using the operator L, this may he written in the form . , 

( l "IT v2 + 2X(1 + v) ] Q 0 = f 6^ } (262)* 

The function Fg is given hy 

Fg = -(1 + V)f^ - §|[l - V 2 + 2 X(l + v)]m^ (263) 

Equation (262) may he compared with the second of equations (g) 
on page 469 of reference 9 • 

Analysis of edge effect for spherical shell.— The special case of 
no distributed surface load and no concentrated load at the apez of the 
shell is obtained hy setting 


F 5 = Fg = 0 


Following again a known procedure from the theory without 
transverse stress deformation, there may he set 


Ql 


\/ sin 0 

_h 

\j sin 0 




( 264 ) 


/ Qi' 


- ^ cot 

\J sin 0 2 


Qi 


\] sin 0 


f » 


\faln 0 


— cot 0 


V 

\/sin 0 


(y cot20 + 1^_SL 

2 ) ,frr 


'Jsin 0 


y (265) 
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with corresponding formulae for 00 * and 00 *'. Introduction of 
equation (2 6k) into equations (257) and (262) gives 


V + (m + \ ~ i cot 2 0^Qi + YTT/ 3 Pl = 0 (266) 


Pl « • - - 1 + 2 cot 2 0 ^ 0 ! - §|[l - V 2 + 2 X(l + v)] Qi = 0 ( 267 ) 

Attention is restricted to the cases for which F^ = Fg = 0 in 

equations (266) and (267) and the problem is considered of the shell 
subject to edge loads (y) Q , 001 ^)o at a section 0 = 00* 

Assuming that cot 0 O is no-t large compared with unityl® and that the 
effect of the edge loads is restricted to a narrow edge zone so 

that | Q]_| « |q 1 "|, | pj « 1 0]_* ' | ) equations (2 66) and (267) maybe 

simplified to 


Ql 


T I _ 


niQl + 


c* 


1 + X /3 


01 = 0 


(268) 



+ 2 X (1 + 



0 


( 269 ) 


Equations (268) and (269) show that the influence of finite E C (X 0 0 ) 
in the edge-effect problem consists, except in extreme circumstances, 
in minor modifications of the results for E c = 00. The quantity 

which represents the influence of finite G c and which is 


~ 2t E f 
^1 ‘ h + t G c 


( 270 ) 


may, however, in practical cases be large compared with unity and not of 
negligible influence on the results. 

■^Ajhen cot » 1 the shell is termed a "shallow" shell which is 
not considered in what follows. 
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Equations (268) and (269) may be compared with equations (l 80 ) 
and (182) for the cylindrical shell. This comparison shows that the 
influence of finite G c in the edge-effect problem is of the same 
nature for the spherical and cylindrical shells. Thus, results of the 
same quantitative nature will he obtainable as in the section on 
cylindrical shells under the headings entitled "Finite circular cyli ndr ical 
shell acted upon by edge moments and forces" and "Semi— inf inite shell acted 
upon by edge bending moment and shear force." 

This work is not herein carried further to specific applications. 

It is apparent that such applications may be worked out with hardly any 
more difficulty than when the effect of the core deformability is not 
taken into account. 


CONCLUDING EEMARKS 


A system of basic equations has been derived for the analysis of 
small-deflection problems of sandwich-type thin shells. This Bystem of 
equations reduces to Love’s theory of thin shells when the transverse 
shear and normal stress deformability of the core of the sandwich is of 
negligible importance. The system of basic equations has been applied 
to a number of specific problems from the theory of plates, circular 
rings, circular cylindrical shells, and spherical shells, and it has 
been found that the effects of both transverse shear and transverse 
normal stress deformation may be of such magnitude that an analysis 
which disregards them gives values for deflections and stresses which 
are appreciably in error. 

Numerical calculations have been in the nature of sample calcu- 
lations, illustrating both the use of the equations and the possible 
effects of using them. Examples have been chosen from the point of 
view of relative simplicity as well as with the thought to illustrate 
most clearly the consequences of the extra deformations which have 
been taken into account. It is unavoidable that, in so doing, some of 
the examples may be of little interest for aircraft structural analysis 
and that some problems may not have been analyzed which would have well 
fitted within the contents of this report and which at the same time 
Wodl'd have been of considerable practical importance. 

The general analysis has been restricted by the following two 
order-of-magnitude relations: (l) t/h « 1 and (2) tE f /hE c » 1 

where t is the face-layer thickness, h is the core-layer thickness, 
E-f is the elastic modulus of the isotropic face— layer material, 
and E c is the elastic modulus in the transverse direction of the 
core-layer material. Therewith it is felt that .very likely nearly 
all situations have been covered in which the effect of transverse 
core flexibility is of significant practical importance. It is 
evident, however, that if desired the theory could be extended so as 
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to include cases where one or "both of these two order-of-magnitude 
relations are not satisfied. The main limitation of the present 
analysis is the omission of all finite-deflection and instability- 
effects. 


Massachusetts Institute of Technology 
Cambridge, Mass., May 26, 19^7 
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